Environmental  Laboratory  erdc/el tr-02-26 


US  Army  Corps 
of  Engineers® 

Engineer  Research  and 
Development  Center 


Installation  Restoration  Research  Program 

Vapor-Phase  Transport  of  Explosives 
Compounds  in  Soils 

R.  Ravikrishna,  Sally  L.  Yost,  Cynthia  B.  Price,  September  2002 

Charolett  A.  Hayes,  Kalliat  T.  Valsaraj, 
and  James  M.  Brannon 


A 


Approved  for  public  release;  distribution  is  unlimited. 


Contents 


Preface . vii 

1 —  Introduction . 1 

Background . 1 

Modeling  of  Vapor-Phase  Transport  of  Explosives . 2 

Objectives . 2 

2 —  Materials  and  Methods . 3 

Soils . 3 

Transport  of  UXO  Chemical  Signatures  from  Soils . 3 

Flux  chamber . 3 

Experimental  design . 4 

Analytical  methods . 4 

Equilibrium  adsorption  testing . 5 

Reaction  rate  testing . 6 

Transport  of  Explosives  Chemicals  from  Buried  Sources . 6 

Mathematical  Model . 7 

3 —  Results  and  Discussion . . . 10 

Model  Simulations . 10 

Case  1 . 10 

Case  2 . 10 

Case  3 . 11 

Case  4 . 11 

Transport  of  UXO  Chemicals  from  Soils . 11 

Experimental  soil/water  adsorption  constants . 11 

Reaction  rates . 12 

Experimental  results  and  model  comparison . 13 

Effect  of  soil  temperature  on  flux  from  soils . 25 

Explosive  Chemical  Signature  Transport  from  Buried  Sources . 28 

Application  of  the  Model  to  Soil  Samples  Collected  Near  Buried 

Landmines . 32 

4 —  Conclusions . 39 

References . 41 


Appendix  A:  Summary  of  Experimental  Data . A1 

Appendix  B:  Notation . B1 

SF  298 

List  of  Figures 


Figure  1 .  Flux  of  2,4-DNT  from  three  soils . 14 

Figure  2.  Flux  of  2,6-DNT  from  three  soils . 1 5 

Figure  3.  Flux  of  1 ,3-DNB  from  three  soils . 1 6 

Figure  4.  Comparison  of  experimental  and  model  fit  for  cases  1 

and  2  for  2,4-DNT  at  20-  and  100-ml/min  airflow  rate 
from  LAAP-C  soil . 1 8 

Figure  5.  Comparison  of  experimental  and  model  fit  for  cases  1 
and  2  for  2,6-DNT  at  20-  and  100-ml/min  airflow  rate 
from  LAAP-C  soil . 19 

Figure  6.  Comparison  of  experimental  and  model  fit  for  cases  1 

and  2  for  1 ,3-DNB  at  20-  and  100-ml/min  airflow  rate 
from  LAAP-C  soil . 20 

Figure  7.  Experimental  data  for  2,4-DNT,  2,6-DNT,  and  1,3-DNB 

for  all  three  soils  for  case  3 . 22 

Figure  8.  Comparison  of  experimental  fluxes  of  2,4-DNT, 

2,6-DNT,  and  1 ,3-DNB  for  case  3  data  with  LAAP-C 

soil  for  20-  and  100-ml/min  airflow  rates . 24 

Figure  9.  Comparison  of  experimental  and  model  fit  for  2,4-DNT, 
2,6-DNT,  and  1,3-DNB  for  case  4  for  20-  and 
1 00-ml/min  airflow  rates,  LAAP-C  soil,  20-percent  soil 
moisture  with  humid  air . 26 

Figure  1 0.  Effect  of  temperature  on  the  flux  of  2,4-DNT  from 

LAAP-C  soil . 27 

Figure  1 1 .  Comparison  of  experimental  and  modeled  soil  concen¬ 
trations  of  2,4-DNT,  2,6-DNT,  and  1,3-DNB  from  a 
buried  source,  soil  moisture  of  5  percent  with  vapor- 
phase  diffusion . 31 

Figure  12.  Comparison  of  experimental  and  modeled  soil  concen¬ 
trations  of  2,4-DNT,  2,6-DNT,  and  1,3-DNB  from  a 
buried  source,  soil  moisture  of  20  percent  with  aqueous- 
phase  diffusion . 33 


Figure  13.  Comparison  of  experimental  and  model  fit  soil  concen¬ 
trations  of  2,4-DNT,  2,6-DNT,  1 ,3-DNB,  and  TNT  from 
buried  TMA-5  type  landmine.  Fort  Leonard  Wood  site, 

41  days  following  burial . 36 

Figure  14.  Comparison  of  experimental  and  model  fit  soil 

concentrations  of  2,4-DNT,  2,6-DNT,  1,3-DNB,  and 

TNT  from  buried  TMA-5  type  landmine.  Fort  Leonard 

Wood  site,  270  days  following  burial . 37 

Figure  Al.  Experimental  soil  concentrations  of  2,4-DNT  as  a 

function  of  position  with  respect  to  source  at  4-cm  depth . A2 

Figure  A2.  Experimental  soil  concentrations  of  2,4-DNT  as  a 
function  of  position  with  respect  to  source  at  1 5-cm 
depth . A3 

Figure  A3.  Experimental  soil  concentrations  of  2,6-DNT  as  a 

function  of  position  with  respect  to  source  at  4-cm  depth . A4 

Figure  A4.  Experimental  soil  concentrations  of  2,6-DNT  as  a 
function  of  position  with  respect  to  source  at  15-cm 
depth . A5 

Figure  A5.  Experimental  soil  concentrations  of  1,3-DNB  as  a 

function  of  position  with  respect  to  source  at  4-cm  depth . A6 

Figure  A6.  Experimental  soil  concentrations  of  1,3-DNB  as  a 
function  of  position  with  respect  to  source  at  1 5-cm 
depth . A7 

Figure  Al.  Experimental  soil  concentrations  of  1,3,5-TNB  as  a 

function  of  position  with  respect  to  source  at  4-cm  depth . A8 

Figure  A8.  Experimental  soil  concentrations  of  1,3,5-TNB  as  a 
function  of  position  with  respect  to  source  at  1 5-cm 
depth . A9 

List  of  Tables 


Table  1.  Physical  Characteristics  of  Soils . 3 

Table  2.  Relevant  Properties  of  Explosives  Compounds . 5 

Table  3.  Initial  Soils  Loading  with  Explosives  Compounds . 5 

Table  4.  Soil  Properties  for  Simulation . 1 0 

Table  5.  Soil- Water  Adsorption  Coefficients  for  2,4-DNT,  2,6- 

DNT,  and  1 ,3-DNB  in  three  soils  at  24  °C . 12 

Table  6.  First-Order  Reaction  Rate  Constants  for  Yokena  Clay 

Soil 


13 


Table  7.  Adjustable  Parameter  for  Comparison  with  Experimental 

Fluxes  from  Cases  1  and  2  for  the  20-ml/min  Airflow 
Rate . 17 

Table  8.  Adjustable  Parameter  for  Comparison  with  Experimental 

Fluxes  from  Cases  1  and  2  for  LAAP-C  Soil  for  the 
1 00-ml/min  Airflow  Rate . 21 

Table  9.  Fit  K$a  and  k\  Values  for  Modeling  Experimental  Data  in 

Case  3  from  LAAP-C  Soil . 25 

Table  10.  Fit  Kaw  Values  to  Model  Experimental  Data  in  Case  4 

from  LAAP-C  Soil . 27 

Table  1 1 .  Adjustable  Parameters  for  Comparison  with  Laboratory 
Flux  Experiments  and  Buried  Chemical  Experiments 
with  LAAP-C  Soil . 32 

Table  12.  Comparison  of  Air  Concentrations  Calculated  from  the 
Buried  Source  Experimental  Data  and  Corresponding 
Model  Simulations . 34 

Table  13.  Model  Fit  Parameters  for  the  Comparison  of  Landmine 

Data  with  Model  Simulations . 38 

Table  A1 .  Initial  and  Final  Concentration  in  Soils  as  a  Function  of 

Temperature  and  Moisture . A10 

Table  A2.  Mass  of  Explosives  from  Buried  Explosives  at  5-Percent 

Moisture  and  4-cm  Burial  Depth . A1 1 

Table  A3.  Mass  of  Explosives  from  Buried  Explosives  at  5-Percent 

Soil  Moisture  and  15-cm  Burial  Depth . A12 

Table  A4.  Mass  of  Explosives  from  Buried  Explosives  at 

20-Percent  Soil  Moisture  and  4-cm  Burial  Depth . A13 

Table  A5.  Mass  of  Explosives  from  Buried  Explosives  at 

20-Percent  Soil  Moisture  and  15-cm  Burial  Depth . A14 


Preface 


The  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC), 
Environmental  Laboratory  (EL),  Vicksburg,  MS,  in  conjunction  with  the 
Gordon  A.  and  Mary  Cain  Department  of  Chemical  Engineering,  Louisiana  State 
University  (LSU),  Baton  Rouge,  LA,  prepared  this  report.  The  experimental  data 
were  generated  at  the  ERDC,  and  the  modeling  was  conducted  at  LSU.  The 
research  was  sponsored  by  Headquarters,  U.S.  Army  Corps  of  Engineers 
(HQUSACE),  Installation  Restoration  Research  Program  (IRRP),  Buried 
Unexploded  Ordnance  Detection  Thrust  Area.  The  Work  Area  was  entitled 
Characterization  and  Mobilization  of  UXO.  Dr.  Clem  Myer  was  the  IRRP 
Coordinator  at  the  Directorate  of  Research  and  Development,  HQUSACE,  and 
the  IRRP  Program  Manager  was  Dr.  M.  John  Cullinane,  EL,  ERDC. 

The  report  was  written  by  R.  Ravikrishna  and  Kalliat  T.  Valsaraj,  LSU,  and 
Sally  L.  Yost,  Cynthia  B.  Price,  Charolett  A.  Hayes,  and  James  M.  Brannon,  EL, 
ERDC. 

This  report  was  reviewed  by  Drs.  Judith  C.  Pennington  and  Patrick  Deliman, 
EL,  ERDC.  The  study  was  conducted  under  the  direct  supervision  of 
Dr.  Richard  E.  Price,  Chief,  Environmental  Processes  and  Engineering  Division, 
EL,  and  Dr.  Edwin  A.  Theriot,  Director,  EL. 

At  the  time  of  publication  of  this  report,  Director  of  ERDC  was  Dr.  James  R. 
Houston.  Commander  and  Executive  Director  was  COL  John  W.  Morris  III,  EN. 

This  report  should  be  cited  as  follows: 

Ravikrishna,  R.,  Yost,  S.  L.,  Price,  C.  B.,  Hayes,  C.  A.,  Valsaraj, 

K.  T.,  and  Brannon,  J.  M.  (2002).  “Vapor-phase  transport  of 
explosives  compounds  in  soils,”  ERDC/EL  TR-02-26, 

U.S.  Army  Engineer  Research  and  Development  Center, 

Vicksburg,  MS. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication ,  or 
promotional  purposes .  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


1  Introduction 


Background 

Military  test  and  training  ranges  contain  large  numbers  of  unexploded 
ordnance  (UXO).  The  detection  and  clearing  of  UXO  at  such  sites  is  currently 
the  Army’s  highest  priority  environmental  restoration  problem.  Chemical  signa¬ 
tures  emanating  from  UXO  potentially  offer  a  means  of  detection,  discrimination 
between  UXO  and  innocuous  clutter,  or  refinement  in  classification  of  detected 
UXO.  The  characterization  and  remediation  of  sites  containing  UXO  using 
currently  available  technology  are  often  hazardous  and  prohibitively  expensive 
due  to  the  extraordinarily  high  level  of  false  detections.  Before  the  utility  of  new 
chemical  detectors  can  be  evaluated,  methods  for  determining  the  strength  of 
explosives  chemical  signatures  under  various  environmental  and  soil  conditions 
are  required. 

The  efficacy  of  chemical  sensors  and  their  potential  usefulness  for  detecting 
buried  UXO  are  difficult  to  determine  without  understanding  how  explosives 
chemical  signatures  are  transported  through  soil.  Soil  conditions  may  result  in  no 
detectable  explosives  signature  at  the  soil  surface  due  to  explosives  degradation 
in  the  soil  (Price,  Brannon,  and  Hayes  1997;  Price  et  al.  2001).  The  probability 
of  detecting  explosives  signatures  is  further  complicated  by  the  low  vapor  pres¬ 
sure  of  2,4,6-trinitrotoluene  (TNT),  hexahydro-l,3,5-trinitro-l,3,5-triazine 
(RDX),  and  octahydro-l,3,5,7-tetranitro-l,3,5,7-tetrazocine  (HMX)  (summarized 
in  McGrath  1995).  The  vapor  pressure  of  RDX  and  HMX  is  so  low  that  signa¬ 
tures  from  spiked  soils  could  not  be  detected  even  with  the  use  of  radioisotopes. 
TNT  forms  a  significant  fraction  of  explosives  in  UXO  because  of  its  use  alone 
or  as  a  component  of  military  formulations  in  munitions.  Manufacturing  impuri¬ 
ties  in  TNT  such  as  2,4-dinitrotoluene  (2,4-DNT),  2,6-dinitrotoluene  (2,6-DNT), 
1,3 -dinitrobenzene  (1,3-DNB),  and  1,3,5-trinitrobenzene  (1,3,5-TNB)  have 
higher  vapor  pressures  than  TNT  and  offer  a  potentially  better  chemical  tracer  for 
detection  of  explosives  from  UXO  (Brannon  and  Pennington,  In  Press;  Phelan 
and  Webb  1997).  The  compounds  2,4-DNT  and  1,3-DNB  were  found  in  higher 
concentrations  in  the  equilibrium  vapor  of  TNT  and  were  transported  more 
rapidly  in  the  soil  (George  et  al.  1999). 

Documentation  exists  on  the  fate  and  transport  of  explosives  residues  from 
contaminated  soils  and  groundwater  under  saturated  conditions  (McGrath  1995; 
Brannon  and  Myers  1997;  Brannon  and  Pennington,  In  Press).  Field  measure¬ 
ments  of  chemical  transport  from  buried  landmines  also  have  been  conducted 
(Jenkins  et  al.  2000;  George  et  al.  1999;  Jenkins,  Leggett,  and  Ranney  1999). 
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However,  descriptions  for  processes  controlling  the  migration  of  explosives 
chemical  signatures  through  soil  from  UXO  and  other  sources  under  unsaturated 
conditions  are  limited  (George  et  al.  1999,  Phelan,  Rodacy,  and  Barnett  2001; 
Phelan  and  Webb  1997)  and  must  be  understood  to  develop  mathematical  models 
to  predict  their  vapor-phase  movement.  These  models  will  aid  in  the  calibration 
and  evaluation  of  existing  and  future  sensors  by  attempting  to  predict  the  strength 
of  chemical  signatures  at  the  detector. 


Modeling  of  Vapor-Phase  Transport  of  Explosives 

The  fate  and  transport  of  explosives  in  the  air-filled  pores  within  soil  affect 
both  the  potential  detection  of  buried  ordnance  by  chemical  sensors  and  vadose 
zone  transport  of  explosives  residues  (McGrath  1995).  The  transport  character¬ 
istics  of  explosives  vapors  through  soils  affect  the  sensitivity  levels  that  sensors 
must  attain  to  detect  chemical  signatures  for  a  given  munitions  depth.  Transport 
characteristics  also  affect  how  chemical  concentrations  will  migrate  in  the  vadose 
zone  independently  of  water-phase  transport.  Experimental  and  modeling 
methods  for  examining  the  flux  of  polyaromatic  hydrocarbons  from  sediment 
into  the  air  have  been  developed  (Valsaraj  et  al.  1999,  1997;  Ravikrishna  et  al. 
1998).  Results  have  shown  that  transport  from  the  sediment  particles  to  the  air 
are  affected  by  moisture  content,  air  relative  humidity,  air  velocity,  and  tempera¬ 
ture.  Documents  have  also  been  published  concerning  the  explosives  source  term 
for  UXO  (Price  et  al.  2000)  and  a  conceptual  model  for  fate  and  transport  of 
explosives  through  soils  (Brannon  et  al.  1999). 

The  modeling  and  experimental  methods  used  for  examining  the  flux  of 
polyaromatic  hydrocarbons  can  be  modified  and  adapted  for  examining  the 
transport  of  explosives  from  soils  under  various  conditions  of  soil  moisture, 
temperature,  and  relative  humidity.  These  parameters  control  the  flux  of  non¬ 
polar  organic  compounds  from  sediment  particles  into  the  air  and  provide  a  basis 
for  evaluation  of  explosives  residues  vapor  transport  (Valsaraj  et  al.  1 999,  1 997; 
Ravikrishna  et  al.  1998). 


Objectives 

The  objectives  of  this  study  were  to  quantify  chemical  signature  transport 
through  soils  under  various  environmental  conditions  in  unsaturated  soils  and  to 
develop  a  model  for  chemical  signature  transport  in  unsaturated  soils. 
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2  Materials  and  Methods 


Soils 

Two  aquifer  soils  obtained  from  Louisiana  Army  Ammunition  Plant, 
Shreveport,  LA  (LAAP-C  and  LAAP-D),  and  one  surface  soil  (Yokena  clay) 
from  the  Mississippi  River  floodplain  were  used  as  test  matrices  in  the  experi¬ 
ments  to  study  UXO  chemical  signature  transport.  These  soils  provided  a  wide 
range  of  physical  and  chemical  characteristics  (Table  1).  Each  soil  was  air-dried, 
ground,  and  sieved  before  storing  at  room  temperature.  After  spiking,  soils  were 
placed  on  a  pan  in  a  thin  layer  under  a  hood  to  allow  the  acetone  to  evaporate, 
then  tumbled  overnight  to  ensure  complete  mixing. 


Table  1 

Physical  Characteristics  of  Soils 

Soil  Property 

LAAP-C 

LAAP-D 

Yokena  Clay 

Percent  sand 

77 

27 

13.75 

Percent  silt 

11 

41 

37.54 

Percent  clay 

12 

32 

48.75 

Textural  classification 

Sandy  loam 

Clay  loam 

Clay 

Percent  total  organic  carbon 

0.08 

0.20 

2.4  | 

Cation  exchange  capacity,  meq/100  g 

6.6 

15.5 

38.9 

Bulk  density  g/cm3 

1.43 

0.88 

0.86 

Transport  of  UXO  Chemical  Signatures  from  Soils 

Flux  chamber 

These  experiments  were  designed  to  measure  and  model  evaporative  fluxes 
of  the  test  explosive  chemicals  from  soils.  The  test  soils  were  spiked  with  the 
chemicals  and  placed  in  a  flux  chamber  with  air  flowing  across  the  surface  of  the 
soil.  A  modified  version  of  a  flux  chamber  designed  by  Spencer  et  al.  (1979) 
was  used.  The  two-part  chamber  was  constructed  of  anodized  aluminum.  The 
bottom  portion  held  soil  at  a  depth  of  2.54  cm  with  a  surface  area  of  30  cm2.  The 
top  portion  was  designed  with  channels  to  provide  a  2-mm  air  space  over  the 
sediment  well  to  allow  uniform  airflow  across  the  soil  surface.  O-rings  and 
threaded  fasteners  were  used  to  seal  the  compartments  together  for  an  airtight  fit. 
Laboratory  house  air,  channeled  through  Gilmont  flow  meters  to  ensure  constant 
flow  rate,  was  passed  through  the  chamber  entrance  ports  over  the  surface  of  the 
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sediment.  Stainless  steel  air  sampling  traps  were  made  from  tubing  and  con¬ 
tained  0.2  g  Tenax  TA  (25-30  mesh).  Traps  were  attached  to  the  exit  ports  of  the 
chambers  to  collect  explosives  compounds  released  into  the  air  stream. 


Experimental  design 

The  experiments  were  designed  to  compare  volatile  emissions  under  varying 
air  relative  humidity  (r.  h.)  (0  and  98  percent),  soil  moisture  content  (5  and 
20  percent),  soil  temperature  (14  and  24  °C)  and  airflow  rate  (20  and 
100  ml/min).  The  100-ml/min  experiments  were  conducted  at  5-  and  20-percent 
soil  moisture  contents  and  with  0-  and  98-percent  r.  h.  conditions,  but  with  one 
soil  only,  LAAP-C.  Air  humidity  was  controlled  using  an  in-line  bubble  trap  to 
add  moisture  vapor  to  the  passing  air.  Moisture  contents  of  each  soil  were 
determined,  and  deionized  water  was  added  to  obtain  the  desired  soil  moisture 
content.  A  recirculating  water  bath  (Remcor  Products  Company,  IL,  model 
CFF550)  was  set  to  maintain  a  constant  14  °C  for  the  low  temperature  experi¬ 
ment,  while  room  temperature  provided  the  desired  24  °C  ambient  temperature. 

Humid  airflow  was  established  at  the  two  desired  flow  rates  of  20  and 
100  ml/min,  and  passed  over  the  soils  for  a  21 -day  period.  The  humid  air  was 
then  switched  to  dry  air  for  another  2 1 -day  period.  Air  samples  were  taken  at  24, 
72,  168,  336  and  504  hr  by  replacing  the  Tenax  trap  with  a  fresh  trap.  The  traps 
were  extracted  with  5  ml  acetonitrile  and  analyzed  for  the  mass  of  contaminant 
trapped  (A m)}  Knowing  the  duration  of  sampling  (A/)  and  the  soil-air  interface 
area  of  the  chamber  (Jc),  the  flux  rates  (NA)  were  determined  from  NA  =  (Am/Ac 
At). 


In  all  flux  measurement  experiments,  the  soils  were  spiked  with  a  mixture  of 
TNT,  2,4-DNT,  2,6-DNT,  1,3,5-TNB,  and  1,3-DNB.  Samples  of  800  g  of  the 
different  soils  were  spiked  with  50  ml  of  an  acetone  solution  containing  approx¬ 
imately  1 0  mg/L  of  each  of  the  explosives.  The  spiked  soils  were  placed  under  a 
hood  to  allow  the  acetone  to  evaporate,  then  tumbled  overnight  to  ensure  com¬ 
plete  mixing.  Table  2  lists  the  physicochemical  properties  of  the  compounds 
considered  in  this  study.  Table  3  lists  the  loading  of  three  of  the  compounds 
spiked  in  the  flux  experiments.  These  three  compounds  were  chosen  for 
comparison  with  the  model  since  their  fluxes  were  consistently  measurable. 


Analytical  methods 

Soils  in  the  flux  chambers  were  analyzed  initially  and  at  the  end  of  each 
experiment  to  determine  if  explosives  degradation  had  occurred.  Air  sampling 
traps  were  extracted  for  24  hr  with  acetonitrile  before  analysis.  At  the  conclusion 
of  an  experiment,  the  plate  at  the  top  of  the  flux  apparatus  was  rinsed  with 
acetonitrile  to  determine  if  any  residues  remained  on  the  chamber  surface. 
Analyses  of  soils  were  performed  by  high-performance  liquid  chromatography 
(HPLC)  using  SW-846  Method  8330  (U.S.  Environmental  Protection  Agency 


1  For  convenience,  symbols  are  listed  in  the  notation  (Appendix  B). 
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Table  2 

Relevant  Properties  ol 

Explosives  Compounds 

Property 

1,3-DNB 

2,4-DNT 

2,6-DNT 

vrrmm 

Molecular  weight 

168.1 

182.1 

182.1 

227.1 

213.1  | 

Diffusivity  in  air,(1)  cm2/s 

0.073 

0.067 

0.067 

0.064 

0.068 

Henry’s  Law  constant,1 [V 
dimensionless 

Boundary  layer  mass  transfer 
coefficient,  ka,  m/s 

20-ml/min  airflow  rate 

1 00-ml/min  airflow  rate 

0.00016 

0.000357 

0.00015 

0.000337 

0.00015 

0.000337 

0.000146 

0.000327 

0.000152 

0.000340 

1  From  summary  in  McGrath  (1995). 

2  From  Thibodeaux  (1996). 

Table  3 

1  Initial  Soils  Loading  with  Explosives  Compounds 

Loading,  mg/kg  fl 

Soil 

1,3-DNB 

2,4-DNT  2,6-DNT  1 

1  20-ml/min  experiments  || 

LAAP-C 

7.4 

10.5 

8.4 

LAAP-D 

9.2 

11.0 

9.0 

Yokena  clay 

9.1 

9.3 

9.4 

1  100-ml/min  experiments  || 

1  LAAP-C 

9.74 

10.4 

9.36  | 

(U.S.  EPA)  1994).  Analyses  of  sampling  traps,  acetonitrile  extracts,  and 
acetonitrile  rinses  of  flux  chamber  surfaces  were  preformed  by  gas  chromato¬ 
graph  using  SW  846,  update  IV-B  Method  8095  (U.S.  EPA  2000). 


Equilibrium  adsorption  testing 

One  of  the  most  critical  parameters  in  the  transport  of  chemicals  in  sub¬ 
surface  soils  is  equilibrium  partitioning  of  the  chemical  between  soil  and  water. 
Equilibrium  adsorption  testing  was  conducted  with  LAAP-C,  LAAP-D,  and 
Yokena  clay  soils  with  a  1:4  ratio  of  soil  to  water  (4  g  of  soil  to  16  ml  water)  and 
five  concentrations  of  contaminant  (10,  7.5,  5,  2.5,  and  1  pg/ml).  The  LAAP-C 
test  was  spiked  with  2,6-DNT;  the  LAAP-D  test  was  spiked  witii  2,6-DNT  and 

2.4- DNT  in  separate  runs;  the  Yokena  clay  was  spiked  with  a  mixture  of 
contaminants  (TNT,  2,4-DNT,  1,3,5-TNB,  1,3-DNB,  and  2,6-DNT).  Tests  were 
placed  on  a  reciprocating  shaker  at  280  excursions  per  minute  for  24  hr.  After 
shaking,  tests  were  centrifuged  at  7,000  relative  centrifugal  speed  for  30  min. 
The  aqueous  phase  was  removed  and  frozen  until  analyzed.  Analytes  included 

2.4- DNT,  2,6-DNT,  1,3-DNB,  1,3,5-TNB,  and  2,4,6-TNT,  2,6-DANT, 

2.4- DANT,  and  4,4-Azoxy.  The  aqueous  samples  from  the  Yokena  clay  were 
preserved  with  ethylene  diamine  tetraacetic  acid  to  a  final  concentration  of 

5  pg/ml.  The  testing  was  carried  out  in  duplicate  for  each  experimental  unit. 
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Reaction  rate  testing 

Disappearance  or  degradation  due  to  reaction  within  the  soil  matrix  is 
another  important  parameter  in  the  subsurface  transport  of  chemicals.  These 
reactions  may  be  due  to  metabolism  by  microbial  organisms  present  in  the  soil  or 
reactions  on  the  soil  surface.  The  three  soils  (Yokena  clay,  LAAP-C,  and 
LAAP-D)  were  spiked  with  a  mixture  of  TNT,  1,3-DNB,  2,4-DNT,  2,6-DNT, 
and  1,3,5-TNB  at  a  concentration  of  1 0  mg/kg.  Spiking  was  done  by  misting  the 
soil  with  the  contaminant  solution  and  then  tumbling  for  24  hr.  After  spiking,  the 
soils  were  divided  into  batches  and  brought  to  two  different  moistures,  5  percent 
and  20  percent  by  weight.  Samples  were  weighed  in  5-g  increments  and  loaded 
into  20-ml  amber  glass  vials.  Sufficient  vials  were  loaded  so  that  each  soil  was 
duplicated  at  four  different  conditions;  5-percent  moisture  at  14  °C,  5-percent 
moisture  at  24  °C,  20-percent  moisture  at  14  °C,  and  20-percent  moisture  at 
24  °C.  The  LAAP-C  and  LAAP-D  treatments  were  sampled  at  7,  14,  and  28  days, 
while  the  Yokena  clay  treatments  were  sampled  at  6  hr,  20  hr,  2  days,  7  days,  and 
14  days.  In  addition  to  the  time-series  sampling,  initial  samples  were  taken  of 
each  soil  at  the  5-  and  20-percent  moisture  level  for  HPLC  analysis. 


Transport  of  Explosives  Chemicals  from  Buried 
Sources 

Another  series  of  experiments  was  conducted  to  evaluate  the  transport  of 
explosives  from  a  controlled  system  similar  to  that  of  buried  munitions.  To 
determine  the  effects  of  soil  moisture  and  burial  depth  on  chemical  signature 
movement,  four  tests  were  conducted  using  soil  at  5-  and  20-percent  moisture 
contents  and  explosives  buried  at  depths  of  15  cm  (6  in.)  and  4  cm  (1 .6  in.)  below 
the  soil  surface.  Soil  obtained  from  LAAP-C  was  used  for  the  four  soil  condi¬ 
tions  tested.  Bulk  soil  samples  were  brought  to  5-  and  20-percent  moisture 
content  with  distilled  deionized  water.  Four  37.85-L  (10-gal)  containers  com¬ 
posed  of  high-density  polyethylene  with  measurements  of  38.5  cm  height  by 
39.4  cm  outside  diameter  (OD)  were  used.  Packing  the  soil  into  the  container 
was  necessary  to  obtain  vertical  cores  that  could  be  sectioned  incrementally' 

The  explosives  were  introduced  by  adding  a  mixture  of  100  mg  each  of  TNT, 
1,3-DNB,  1,4,5-TNB,  2,4-DNT,  and  2,6-DNT  to  the  surface  of  soil  15  cm  deep 
through  a  stainless  steel  tube  placed  centrally  in  each  container.  Soil  in  one 
container  was  at  5-percent  moisture;  soil  in  the  other  was  at  20-percent  moisture. 
Approximately  2.54  cm  (1  in.)  of  soil  was  packed  around  the  tube.  The  tube  was 
then  removed  and  the  explosives  were  visually  inspected  to  ensure  that  they  had 
remained  in  place.  Soil  was  then  packed  over  and  around  the  explosive  plug. 
Additional  soil  was  added  and  packed  to  a  depth  of  1 5  cm  over  the  mixture 
giving  a  total  depth  of  30  cm.  Identical  procedures  were  used  to  place  the  explo¬ 
sives  in  the  remaining  two  containers  with  the  exception  that  the  soils  were  ini¬ 
tially  packed  to  a  depth  of  26  cm.  An  additional  4  cm  of  soil  was  packed  over 
the  compounds. 

The  containers  were  sealed  and  allowed  to  sit  at  room  temperature  (-24  °C) 
for  120  days.  Thin-walled  stainless  steel  tubing,  1.9  cm  (0.75  in.)  OD,  was  used 
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to  sample  the  soils.  Eight  horizontal  cores  both  directly  over  or  to  the  side  of  the 
explosives  mixture  were  made  in  the  20-percent  moisture  tests.  The  stainless 
cores  were  sectioned  every  1.5  cm  from  the  soil  surface  giving  approximately 
19  vertical  samples  per  horizontal  core.  Due  to  the  looser  consistency  of  the 
5-percent  moisture  soil,  only  six  horizontal  cores  were  taken  in  each  container  for 
approximately  17  vertical  segments  for  each  core.  Ten  milliliters  of  acetonitrile 
was  added  to  each  section,  and  samples  were  stored  in  the  freezer.  Twenty 
samples  from  each  container  were  selected  and  analyzed  using  SW-846  8095 
(U.S.  EPA  2000). 

Mathematical  Model 

A  general  model  describing  the  contaminant  fate  and  transport  in  porous 
media  is  a  diffusion-reaction  equation: 

rlC1 

Rf~f-  =  ^*(D^CA)-(~rA)  (1) 


where 


Rf  =  retardation  factor 

CA  =  concentration  of  chemical  in  pore  air 


t  =  time 


Deff  =  effective  diffusivity,  m2/s 
-rA  -  reaction  term 

Diffusion  in  the  soil  pore  air  is  assumed  to  be  the  dominant  mechanism  of  mass 
transport  within  the  soil.  The  reaction  term,  -rA,  is  added  for  generality  to  accom¬ 
modate  disappearance  of  chemical  that  may  occur  due  to  biodegradation  or 
reaction  within  the  soil.  The  model  also  neglects  chemical  loss  by  the  mecha¬ 
nism  of  “wicking”  or  vapor-phase  advective  transport  of  the  chemical  that  may 
occur  as  a  result  of  water  evaporation. 

If  only  the  vertical  transport  in  the  z-direction  is  considered,  Equation  1  is 
simplified  to  Equation  2,  which  is  the  governing  differential  equation  for  the  pore 
air  concentration  of  the  contaminant,  CA  ( z,t ),  as  a  function  of  depth,  z,  and  time, 
t.  Here  it  is  assumed  that  the  rate  of  degradation  is  first-order  and,  hence,  -rA  = 
k\CA. 


^C±  = 

dt 


A#  d2C’A 
Rf  dz2 


kxCA 


(2) 
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where  k\  is  the  first-order  constant  for  disappearance  rate  of  chemical  in  soil. 
The  following  boundary  conditions  are  applicable  when  a  contaminated  soil  is 
exposed  to  air  as  applied  to  the  experimental  data  from  the  flux  experiments 
conducted  with  surface  soils.  At  the  soil-air  interface,  a  thin  film  offers  a  mass 
transfer  resistance.  The  boundary  condition  here  at  z  =  0  (surface)  is  described 
by: 


k  Q  =d  ^£a 
aA  dz 


(3) 


where  ka  is  the  air-side  mass  transfer  coefficient.  At  depths  far  away  from  the 
soil-air  interface,  z  =  oo,  the  chemical  concentration  in  the  soil  is  relatively 
unaffected  by  the  transport  dynamics  at  the  air-soil  interface;  therefore,  this 
semiinfinite  boundary  condition  can  be  set  as  CA  =  CA  where  C°A  is  the  initial 
equilibrium  concentration  of  chemical  in  pore  air. 

The  analytical  solution  to  Equation  2  (Choy  and  Reible  2000)  gives  CA(z,t). 
The  flux  is  calculated  by  multiplying  the  surface  overall  mass  transfer  coeffi¬ 
cient,  ka,  and  the  vapor-phase  contaminant  concentration  at  the  surface,  CA(0,t), 
and  the  resultant  expression  for  flux,  NA  (/),  is  shown  in  Equation  4. 


N A(t)  =  C°Aka e~k'’  ■  exp 


JlL 

KD«RfJ 


erfc 


(4) 


where  erfc  is  the  complementary  error  function.  In  Equation  4,  CA ,  the  initial 
equilibrium  concentration  of  the  contaminant  in  the  pore  air,  is  given  by 
fV°  Pj  /  Rf  where  W°  is  the  initial  loading  of  the  contaminant  (kg/kg)  and  p*  is 
the  bulk  density  of  the  soil  (kg/m3). 

The  model  described  by  Equation  1  assumes  local  equilibrium  with  respect  to 
the  contaminant  in  the  pore  spaces.  Further,  the  surface  flux  depends  on  the  soil- 
side  and  air-side  resistances  offered  to  transport  of  the  contaminant  (Valsaraj  et 
al.  1999).  The  soil-side  resistance  is  dependent  on  the  retardation  of  the  contam¬ 
inant  on  the  soil  surface  due  to  adsorption  and  the  effective  diffusivity  of  the 
contaminant  in  the  porous  media.  The  retardation  factor  is  given  by  ya  +  p  fCSA, 
where  K$A  is  the  soil-air  partition  constant  (L/kg)  and  ya  is  the  air-filled  porosity 
(m3/m3).  The  retardation  factor  is  proportional  to  the  partition  constant  of  the 
contaminant  between  soil  and  air  (KSA).  In  the  absence  of  direct  measurements, 
KsA  can  be  estimated  for  wet  soils  (>5-percent  soil  moisture  content)  using 
Equation  5  (Valsaraj  et  al.  1 997): 


K 


SA 


kaw 


(5) 
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where 


Ksw=  soil-water  partition  constant,  L/kg 

Kaw=  air-water  partition  constant,  Henry’s  Law  constant,  molar 
concentration  ratio,  dimensionless 

The  soil-water  partition  constants,  Ksw,  were  directly  measured  for  the  com¬ 
pounds  of  interest  in  separate  batch  experiments.  In  the  absence  of  measured 
data,  the  values  of  the  Henry’s  Law  constant  for  the  compounds  of  interest  were 
obtained  from  the  ratio  of  the  saturated  vapor  pressure  and  solubility  of  explo¬ 
sives  in  water.  The  uncertainty  in  the  Henry’s  Law  constant  values  directly 
translates  to  the  uncertainty  in  the  soil-air  partition  constant. 

The  effective  diffusivity,  Deg,  was  computed  using  the  expression,  DA  yo10/3  / 
yT2,  where  DA  is  the  diffusivity  of  contaminant  A  in  air,  ya  is  the  air-filled 
porosity,  and  yj-is  the  total  porosity  (Choy  and  Reible  2000).  The  values  of  air- 
filled  porosity  used  in  the  simulation  were  calculated  using  the  measured  values 
of  the  initial  soil  moisture  content,  total  soil  porosity,  and  soil  bulk  density.  The 
mass  transfer  coefficient,  ka,  quantifies  the  air-side  resistance  offered  by  the  film 
for  mass  transfer  between  the  soil  surface  and  the  air  flowing  over  the  surface. 
The  ka  was  estimated  using  the  boundary  layer  theory,  0.664  Re0  33  Sc° 5  ( DAld ), 
where  Re  is  the  Reynolds  Number,  Sc  is  the  Schmidt  number,  and  d  is  the  length 
of  the  airflow  path.  The  Re  is  given  by  d-v/v,  where  v  is  the  average  flow 
velocity,  v  is  the  kinematic  viscosity  of  air,  and  Sc  is  given  by  v/DA.  All 
parameters  were  assumed  to  be  constant  during  the  simulation. 

Wicking  can  be  an  important  transport  mechanism  when  soil  moisture  is  con¬ 
tinuously  replenished  from  beneath  the  top  layers  of  soil  during  surface  drying. 
The  laboratory  experiments  described  in  this  manuscript  were  performed  with  a 
very  thin  layer  of  soil  (2.54  cm)  that  had  low  moisture  content  (5  or  20  percent) 
and  no  water  replenishment  from  below. 

The  losses  from  any  possible  wicking  can  be  estimated  by  kaCh  •  Cair°,  where 
kadv  is  the  advective  velocity  that  may  be  induced  by  wicking,  and  Ca,r°  is  the 
equilibrium  vapor-phase  concentration  of  the  chemical  in  the  soil.  The  advective 
velocity  of  pore  air  saturated  with  water  vapor  was  calculated  from  a  previous 
study  (Valsaraj  et  al.  1997)  and  corrected  for  the  difference  in  airflow  rates  used 
in  both  studies  based  on  soil  to  air  mass  transfer  coefficients  in  the  literature 
(Thibodeaux  1996),  which  was  0.034  cm3/cm2/hr.  The  advective  fluxes  thus 
computed  for  the  three  chemicals  discussed  in  this  report  ranged  from  7.4  x  10‘6 
to  0.025  ng/cm2/hr.  The  experimental  fluxes  measured  at  times  below  72  hr 
(time  period  when  steady-state  water  evaporation  is  expected  to  occur)  ranged 
from  0.0083  to  17.2  ng/cm2/hr.  The  ratios  of  the  computed  advective  fluxes  to 
the  corresponding  measured  experimental  fluxes  ranged  from  0.00044  to  0.05.  In 
the  worst  case,  the  advective  fluxes  were  5  percent  of  the  measured  fluxes.  This 
analysis  shows  that  for  chemicals  that  exhibit  very  low  vapor  pressures  and  very 
high  partition  constants,  wicking  was  not  a  significant  mechanism  for  vapor- 
phase  transport  and  was,  therefore,  neglected  in  the  transport  model  described 
herein. 
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3  Results  and  Discussion 


Model  Simulations 

Experimental  data  analyzed  using  the  model  was  classified  under  four  cases 
(Ravikrishna  et  al.  2002)  on  the  basis  of  the  experimental  conditions  (initial  soil 
moisture  content  and  relative  humidity  of  air  passing  over  the  soil  surface).  The 
porosity  values  used  in  each  case  of  the  simulation  are  summarized  in  Table  4. 


Table  4 

Soil  Propert 

ies  for  Simulation 

Experiment 

Initial  Soil 
Moisture,  % 

Air  Relative 
Humidity 

Air-Filled 
Porosity,  Ya 

Water-Filled 
Porosity,  yw 

Total 

Porosity,  yt 

Case  1 

5 

humid 

0.23 

0.08 

0.31 

Case  2 

20 

dry 

0.16 

0.15 

0.31 

Case  3 

5 

dry 

0.31 

0 

0.31 

Case  4 

20 

humid 

0 

0.31 

0.31 

Case  1 

The  initial  moisture  content  of  the  soil  was  5  percent  and  the  relative 
humidity  of  the  air  passing  over  the  soil  surface  was  1 00  percent.  With  humid  air 
passing  over  the  soil,  no  evaporation  was  expected.  The  soil  pore  air  was 
expected  to  retain  the  initial  moisture  content.  The  initial  air-filled  porosity  was 
calculated  to  be  0.23  and  was  maintained  at  that  value  throughout  the  simulation 
period.  At  5-percent  soil  moisture  content,  the  soil  was  considered  wet  and, 
therefore,  the  KSA  was  estimated  from  Ksw  and  KAW  as  per  Equation  4. 


Case  2 

The  initial  moisture  content  of  the  soil  was  20  percent,  and  the  relative 
humidity  of  the  air  passing  over  the  soil  surface  was  0  percent.  Since  the  air 
passing  over  the  soil  was  dry,  evaporation  was  expected.  Soil  moisture  filled 
almost  all  of  the  soil  pores  initially,  but  with  moisture  loss  to  air,  the  air-filled 
porosity  was  expected  to  increase.  Earlier  reports  of  evaporation  rates  from 
similar  experiments  (Valsaraj  et  al.  1997)  showed  that  initially  water  loss  occurs 
very  rapidly  and  complete  water  loss  takes  much  longer.  The  air-filled  porosity 
would  increase  as  a  function  of  time  during  this  period,  but  in  the  absence  of  any 
transient,  measured  data,  the  air-filled  porosity  was  set  as  an  average  of  the  initial 
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(zero)  and  the  final  expected  value  (total  porosity).  For  the  duration  of  the 
experiment,  the  soil  was  presumed  to  retain  enough  moisture  to  be  considered 
wet,  and  the  Ksa  was  estimated  from  Ksw  and  KAW  as  per  Equation  4. 


Case  3 

The  initial  moisture  content  of  the  soil  was  5  percent,  and  the  relative  humid¬ 
ity  of  the  air  passing  over  the  soil  surface  was  0  percent.  The  moisture  from  the 
pore  air  space  was  expected  to  decrease  within  a  very  short  period  of  time  and, 
therefore,  the  value  of  the  air-filled  porosity  was  set  equal  to  that  of  the  total 
porosity  throughout  the  simulation  period.  The  rapid  drying  creates  the  possi¬ 
bility  of  dry-off  of  the  soil  surface.  In  this  case,  Ksa  cannot  be  approximated  as 
Ksw! He  (where  He  =  Henry’s  Law  constant)  and  was  not  directly  measured 
either.  Ksa  was  used  as  an  adjustable  parameter  in  the  model  for  this  case. 


Case  4 

The  initial  moisture  content  of  the  soil  was  20  percent,  and  humid  air  was 
passed  over  the  soil  during  the  entire  course  of  the  experiment.  Initially,  the  soil 
pore  spaces  were  all  filled  with  water,  and  this  condition  persisted  throughout  the 
entire  experiment.  This  implies  that  the  diffusion  through  pore  water  is  the 
predominant  transport  mechanism  rather  than  the  diffusion  through  pore  air. 
Equation  2  remains  the  governing  differential  equation,  but  CA  is  the  pore  air 
concentration  instead  of  pore  water  concentration.  The  initial  pore  water  concen¬ 
tration,  C°A ,  is  given  by  I  Ksw,  where  Ksw  is  the  measured  soil-water 

equilibrium  partition  constant.  At  the  soil-air  interface,  a  correction  factor  of  KAW 
(the  air-water  equilibrium  constant,  or  the  Henry’s  Law  constant)  is  applied  to 
the  overall  mass  transfer  coefficient,  ka,  to  account  for  the  evaporation  at  the  soil- 
air  interface.  The  effective  diffusivity,  DCff,  is  recomputed  using  water-filled 
porosity  and  the  diffusivity  of  the  chemical  in  water.  The  retardation  factor,  Rf, 
was  computed  as  yw  +  Ksw  Pa,  where  yw  =  water-filled  porosity. 


Transport  of  UXO  Chemicals  from  Soils 

Experimental  soil/water  adsorption  constants 

The  two  aquifer  soils  from  LAAP  were  high  in  sand,  ranging  from  27  to 
77  percent  sand  and  low  in  organic  carbon  (Table  1).  Silt  and  clay  were  present 
in  all  samples,  although  in  lower  amounts  in  the  LAAP-C  soil.  Cation  exchange 
capacity  (CEC)  was  also  relatively  low,  ranging  from  6.6  to  15.5  meq/100  g.  The 
Yokena  clay  surface  soil  was  high  in  clay,  organic  carbon,  and  CEC  compared  to 
the  aquifer  soils. 

The  experimentally  measured  soil/water  partition  constants  (Ksw)  for 
2,4-DNT,  2,6-DNT,  and  1,3-DNB  are  listed  in  Table  5.  Soil/water  partition 
constants  for  TNT  are  presented  elsewhere  (Pennington  et  al.  1999).  The  KSw 
values  for  the  Yokena  clay  surface  soil  were  larger  than  those  for  the 
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Table  5 

Soil-Water  Adsorption  Coefficients  ( Ksw ,  L/kg)  for  2,4-DNT,  2,6-DNT, 
and  1,3-DNB  in  Three  Soils  at  24  °C 

Soil 

Compound 

Ksw 

R2 

LAAP-C  (aquifer  soil) 

2,4-DNT1 

0.67 

0.85  | 

LAAP-D  (aquifer  soil) 

2,4-DNT 

1.67 

0.75 

Yokena  clay 

2,4-DNT 

12.5 

0.95 

LAAP-C  (aquifer  soil) 

2,6-DNT 

0.96 

0.96 

LAAP-D  (aquifer  soil) 

2,6-DNT 

1.83 

0.88 

Yokena  clay 

2,6-DNT 

5.96 

0.99 

LAAP-C  (aquifer  soil) 

1,3-DNB1 

0.32 

0.59 

LAAP-D  (aquifer  soil) 

1,3-DNB1 

No  significant  adsorption 

Yokena  clay 

1,3-DNB 

17.7 

0.95 

LAAP-C  (aquifer  soil) 

1,3,5-TNB1 

0.49 

0.99 

LAAP-D  (aquifer  soil) 

1,3,5-TNB1 

0.27 

0.88 

Yokena  clay 

1,3,5-TNB 

Steady-state  concentrations  not  attained 

1 1  Data  from  Pennington  et  al.  (1999).  || 

aquifer  soils.  Previous  work  on  a  common  UXO  fill  explosive,  TNT,  reported  a 
strong  correlation  between  KSw  and  the  CEC  and  clay  content  of  the  soil 
(Pennington  and  Patrick  1990;  Haderlein,  Weissmahr,  and  Schwarzenbach  1996). 
In  the  present  work,  the  KSw  for  all  four  compounds  increased  with  CEC,  clay 
content,  and  organic  carbon  fraction  of  the  soils.  Sorption  of  these  compounds  is 
clearly  higher  in  the  surface  soil  (Yokena  clay)  than  in  aquifer  soils  (LAAP-C 
and  LAAP-D). 

Reaction  rates 

Concentrations  of  explosives  remained  essentially  unchanged  in  the  aquifer 
soils  regardless  of  the  temperature  or  moisture  regime  (Appendix  A,  Tables  A1 
and  A2).  This  held  true  for  1,3,5-TNB,  1,3-DNB,  TNT,  2,4-DNT,  and  2,6-DNT 
in  the  aquifer  soils.  Initial  and  final  explosives  concentrations  for  the  Yokena 
clay  soil  are  summarized  in  Table  A3.  Results  ranged  from  no  change  in  soil 
concentrations  over  the  1 4  days  of  incubation  to  complete  disappearance  of  the 
tested  compound. 

The  rate  of  processes  that  remove  TNT  and  its  transformation  products  from 
solution  was  modeled  using  pseudo-first-order  kinetics,  which  take  the  form 
ln(Co/C)  =  k\»t ,  where  C  is  the  chemical  concentration  of  the  reacting  substance, 
(mg/L),  Co  is  the  concentration  of  the  reactant  at  time  0,  k\  is  the  pseudo-first- 
order  rate  constant  (hr1),  and  t  is  time  (hr).  Half-life  (/]/2)  was  calculated  using 
t\i2  =  0.693/Aj. 

The  tested  compounds  generally  disappeared  from  the  soil  under  both  the  5- 
and  20-percent  moisture  regime  at  24  °C  (Table  6)  in  the  Yokena  clay  soil.  Half- 
lives  were  much  shorter  at  higher  soil  moisture  levels.  At  the  5-percent  moisture 
level,  half-lives  of  1,3-DNB  and  2,6-DNT  could  not  be  determined  because  no 
significant  change  in  soil  concentrations  were  observed  over  the  course  of  the 
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Table  6 

First-Order  Reaction  Rate  Constants  (K,  hr'1)  for  Yokena  Clay  Soil 


Compound 

14  °C 

24  °C  ~ 1 

5-Percent 

Moisture 

20-Percent 

Moisture 

5-Percent 

Moisture 

20-Percent  I 

Moisture  | 

K,  hr'1 

p 

K,  hr'1 

p 

K,  hr'1 

p 

K,  hr'1 

p 

1,3,5-TNB 

0.0017 

0.51 

0.027 

0.988 

0.0005 

0.79 

0.093 

0.94 

1,3-DNB 

NSC1 

0.012 

0.91 

NSC 

0.137 

0.92 

TNT 

0.0009 

0.56 

0.033 

0.989 

0.0007 

0.81 

0.122 

0.94 

2,6-DNT 

NSC 

NSC 

NSC 

0.011 

0.99 

2,4-DNT 

NSC 

0.005 

0.99 

0.0005 

0.86 

0.011 

0.91 

1  NSC  =  No  significant  change. 


experiment.  Where  significant  changes  occurred  in  soil  concentrations,  removal 
of  the  explosives  contaminants  generally  conformed  to  first-order  kinetics 
(Table  6).  This  suggests  that  when  modeling  stability  of  these  compounds  in 
soils,  first-order  kinetics  can  be  used  to  describe  removal  that  results  from  the 
sum  of  all  soil  reactions.  Values  of  r2  were  especially  high  under  the  20-percent 
moisture  regime,  exceeding  r  =  0.91  for  compounds  where  soil  concentrations 
decreased.  At  the  20-percent  moisture  regime,  the  tested  compounds  were  more 
stable  at  14  °C  than  at  24  °C.  Half-lives  at  24  °C  were  less  than  half  of  that 
measured  at  14  °C.  These  results  show  that  both  temperature  and  moisture  affect 
the  stability  of  1,3,5-TNB,  1,3-DNB,  TNT,  2,4-DNT,  and  2,6-DNT  in  the 
Yokena  clay  soil.  Higher  temperature  and  moisture  resulted  in  decreased 
stability  of  these  compounds,  which  can  affect  their  signature  transport. 


Experimental  results  and  model  comparison 

Cases  1  and  2, 20-ml/min  flow  rate  experiments.  Figures  1-3  show  the 
comparison  of  experimental  fluxes  and  the  corresponding  model  simulations  for 

2.4- DNT,  2,6-DNT,  and  1,3-DNB,  respectively,  at  24  °C  in  three  different  soils 
for  cases  1  and  2,  and  at  the  airflow  rate  of  20  ml/min.  These  three  compounds 
were  chosen  for  comparison  since  the  data  set  was  most  complete  for  these. 
Figure  1  shows  the  flux  of  2,4-DNT  from  LAAP-C,  LAAP-D,  and  Yokena  clay 
soils.  The  experimental  flux  values  are  represented  by  discrete  symbols  while 
lines  represent  the  simulation  curves.  The  experimental  data  show  the  flux  of 

2.4- DNT  decreasing  gradually  with  time  and  quickly  reaching  steady-state. 
Initially,  the  flux  is  air-phase-resistance  controlled.  Very  quickly  the  flux 
becomes  soil-side-resistance  controlled,  and  diffusion  through  the  soil  pores 
dominates.  Thus,  the  initial  flux  can  be  given  by  NA(t)  =  kaCA  and  as  t  becomes 
large,  the  long-term  flux  is  proportional  to  \ltm.  Any  deviation  from  this 
behavior  is  indicative  of  processes  other  than  diffusion  also  being  significant. 
Figures  2  and  3  show  similar  trends  in  the  experimental  fluxes  of  2,6-DNT  and 
1,3-DNB,  respectively,  for  cases  1  and  2  for  the  three  soils. 

The  model  curves  in  each  case  shown  in  Figures  1-3  were  not  a  priori  simu¬ 
lations,  but  best  fits  of  the  experimental  data  using  one  variable  fit  parameter  - 
the  soil-air  partition  constant,  Ksa ■  The  best  fit  for  each  data  set  was  obtained  by 
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Flux,  ng/cmz/hr  „  Flux,  ng/cm2/hr  Flux,  ng/cm2/hr 


a)  LAAP-Csoil 


b)  LAAP-Dsoil 


Figure  1.  Flux  of  2,4-DNT  from  three  soils.  Comparison  of  experimental  data  and  model  fits  for  cases  1 
and  2  at  20-ml/min  airflow  rate 
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Figure  2.  Flux  of  2,6-DNT  from  three  soils.  Comparison  of  experimental  data  and  model  fits  for  cases  1 
and  2  at  20-ml/min  airflow  rate 
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Flux,  ng/cnF/hr  Flux,  ng/om=/hr  Fta  ng/cm^r 


a)  LAAP-Csoil 


b)  LAAP-D  soil 


c)  Yokena  clay  soil 


0  100  200  300  400  500  600  700 


Time,  hr 

Figure  3.  Flux  of  1 ,3-DNB  from  three  soils.  Comparison  of  experimental  data  and  model  fits  for  cases  1 
and  2  at  20-ml/min  airflow  rate 
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adjusting  the  Ksa  value  until  the  sum  of  squares  quantity  between  the  experi¬ 
mental  data  and  model  fit  was  a  minimum.  The  model  fit  curves  show  a  trend 
similar  to  that  of  the  experimental  data.  The  effectiveness  of  the  model  to 
simulate  the  experimental  fluxes  can  be  evaluated  by  comparing  the  experi¬ 
mentally  determined  and  fit  model  parameters,  in  this  case,  Ksa ■  Table  7  shows 
this  comparison  of  the  KSa  values  for  cases  1  and  2  for  the  three  compounds  and 
the  three  soil  types.  The  average  difference  between  the  fit  and  the  a  priori 
estimation  of  Ksa  is  within  an  order  of  magnitude  for  all  compounds.  The  Ksw 
values  used  to  estimate  Ksa  were  those  measured  in  the  batch  experiments,  and, 
hence,  the  variation  in  Ksa  is  a  direct  result  of  the  variation  in  reported  KAW  values 
(Henry’s  Law  constant). 


Table  7  1 

Adjustable  Parameter  for  Comparison  with  Experimental  Fluxes 
from  Cases  1  and  2  for  the  20-ml/min  Airflow  Rate1 

Partition 

Constant 

LAAP-C 

LAAP-D 

Yokena  Clay 

Case  1  |  Case  2 

Case  1  |  Case  2 

Case  1  |  Case  2 

2,4-DNT  || 

Ksa,  L/kg 
estimated 

8.9e+4 

8.9e+4 

2.2e+5 

2.2e+5 

1.7e+6 

1.7e+6 

Ksa,  L/kg 
fit 

1.2e+5 

1.1  e+5 

5.5e+6 

4.8e+5 

4.2e+6 

4.2e+7 

2,6-DNT  || 

Ksa,  L/kg 
estimated 

1.3e+5 

1.3e+5 

2.4e+5 

2.4e+5 

7.9e+5 

7.9e+5 

Ksa,  L/kg 
fit 

6.4e+3 

2.5e+4 

9.6e+4 

4.7e+4 

3.2e+5 

5.5e+6 

1  1,3-DNB  1 

Ksa,  L/kg 
estimated 

9.9e+3 

9.9e+3 

UK 

NA 

5.5e+5 

5.5e+5 

Ksa,  L/kg 
fit 

6.1  e+4 

6.9e+4 

2.3e+6 

1  NA  =  Experimental  KSw  not  available  and  hence  KSa  not  estimated. 

IMP  =  No  model  fit  performed  since  KSw  measurement  was  not  available. 

ND  =  No  experimental  data. 

Cases  1  and  2, 100-ml/min  airflow  experiments.  Figures  4-6  show  the 
experimental  fluxes  for  2,4-DNT,  2,6-DNT,  and  1,3-DNB,  respectively,  from 
LAAP-C  soil  at  100-  and  20-ml/min  airflow  rates  for  cases  1  and  2.  The  figures 
also  show  a  comparison  with  the  model.  The  model  simulations  were  fit  to  the 
experimental  data  by  varying  the  value  of  KSa .  In  general,  the  experimental 
fluxes  were  higher  in  the  1 00-ml/min  airflow  rate  experiment  than  in  the 
20-ml/min  experiment  for  the  three  compounds.  This  is  expected  since  the 
increased  airflow  rate  implies  a  higher  mass  transfer  coefficient  at  the  soil-air 
interface  that  reduces  the  overall  resistance  for  mass  transfer  from  the  bulk  to  air. 

For  the  last  two  data  points  of  case  2  in  Figures  4-6,  the  experimental  flux  for 
higher  airflow  rate  is  smaller  than  that  for  the  lower  airflow  rate.  This  anomalous 
trend  may  be  explained  by  the  fact  that  in  case  2,  dry  air  is  passing  over  the  sur¬ 
face  of  20-percent  moist  soil.  Over  400  hr,  this  surface  tends  to  dry,  increasing 
the  partition  constant,  which  results  in  lower  flux.  Case  1  represents  constant 
conditions  with  respect  to  the  effective  diffusivity  and  the  partition  constant; 
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Figure  4.  Comparison  of  experimental  and  model  fit  for  cases  1  and  2  for  2,4-DNT  at  20-  and  100-ml/min  airflow  rate  from  LAAP-C  soil 
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Figure  5.  Comparison  of  experimental  and  model  fit  for  cases  1  and  2  for  2,6-DNT  at  20-  and  100-ml/min  airflow  rate  from  LAAP-C  soil 
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Figure  6.  Comparison  of  experimental  and  model  fit  for  cases  1  and  2  for  1 ,3-DNB  at  20-  and  1 00-ml/min  airflow  rate  from  LAAP-C  soil 


therefore,  such  anomalous  behavior  is  not  observed.  The  model  simulations  were 
fit  to  the  experimental  data  using  KSA  as  the  adjustable  parameter.  The  com¬ 
parison  of  the  estimated  and  the  fit  Ksa  values  for  both  flow  rates  is  shown  in 
Table  8.  The  differences  between  the  fit  Ksa  values  in  both  cases  and  the  esti¬ 
mated  Ksa  values  were  discussed  earlier  (Valsaraj  et  al.  1999)  and  were  attributed 
to  the  uncertainty  in  the  reported  values  of  KAW.  All  parameters  used  to  simulate 
the  two  different  fit  cases  were  the  same  except  for  the  air-side  mass  transfer 
coefficient  and  the  fit  Ksa  value.  The  difference  in  the  fit  Ksa  values  between  the 
two  flow  rate  experiments  is  less  than  one  order  of  magnitude.  This  difference 
represents  the  degree  of  validity  of  the  initial  approach  to  fit  the  experimental 
values  by  varying  Ksa-  If  this  approach  is  correct,  the  difference  between  the  fit 
Ksa  values  in  both  cases  (20-  and  1 00-ml/min  airflow  rates)  should  be  insignifi¬ 
cant.  However,  a  small  difference  in  the  fit  Ksa  values  was  observed  for  both 
cases.  The  difference  was  significant  for  case  1,  where  no  change  in  operating 
conditions  over  the  duration  of  the  experiment  occurred. 


Table  8 

Adjustable  Parameter  for  Comparison  with  Experimental  Fluxes, 

KSa  (L/kg),  from  Cases  1  and  2  for  LAAP-C  Soil  for  the  100-ml/min 
Airflow  Rate 

|  Parameter 

Case  1 

Case  2 

j  2,4-DNT 

Estimated 

8.9e+4 

8.9e+4 

Model  fit  -  20  ml/min 

1.2e+5 

9.7e+4 

Model  fit  - 100  ml/min 

1.6e+5 

2.9e+4 

S  2,6-DNT  || 

Estimated 

1 .3e+5 

1.3e+5 

Model  fit  -  20  ml/min 

6.4e+3 

2.5e+4 

Model  fit -100  ml/min 

1.5e+4 

3.2e+4 

1  1,3-DNB  | 

Estimated 

9.9e+3 

9.9e+3 

Model  fit  -  20  ml/min 

6.1e+4 

6.9e+4 

Model  fit  - 100  ml/min 

1.1  e+5 

2.0e+4 

Case  3, 20-  and  100-ml/min  airflow  rate  experiments.  Figure  7  shows  the 
experimental  data  for  case  3  with  the  airflow  rate  at  20  ml/min.  The  fluxes  show 
a  distinct  sharp  downward  gradient  with  time  not  characteristic  of  a  l/tm 
dependence  for  purely  diffusive  transport  as  observed  in  the  cases  1  and  2 
(shown  in  Figures  1-3).  This  behavior  is  uniformly  observed  for  all  three  com¬ 
pounds  and  for  the  three  soil  types.  This  trend  is  not  indicative  of  an  equilibrium 
state  in  the  soil  as  seen  in  Figures  1-3.  It  is  more  representative  of  a  dynamic 
rate-controlled  (kinetic)  phenomenon.  Two  hypotheses  arise  out  of  these 
observations. 

The  first  hypothesis  is  that  the  soil  is  drying  at  a  certain  rate  that  is  dependent 
on  the  airflow  rate.  The  initial  moisture  content  (5  percent)  was  near  the  lower 
limit  of  the  “wet”  soil  criteria  (>5  percent  moisture  content)  for  soil-air  partition 
constants  described  earlier.  The  soil  partition  constant  increases  as  soil  moisture 
decreases  below  5  percent  and  continues  to  decrease  until  the  soil  moisture  level 
reaches  the  “dry”  state  and  is  constant  for  soil  moisture  levels  below  0.1  percent. 
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The  rate  of  decrease  of  the  partition  constant  in  the  damp  zone  between  the  “diy” 
and  the  “wet”  zone  is  not  well  characterized  and  depends  on  the  drying  rate  of 
the  soil,  which  is  dependent  on  the  flow  rate  of  air  and  initial  moisture  content. 
Therefore,  decreasing  flux  measurements  might  represent  successive  stages  of 
increasing  partition  constant  as  the  soil  is  still  undergoing  loss  of  moisture  in  the 
“damp”  regime  due  to  the  low  moisture  content  and  the  low  airflow  rate 
(20  ml/min).  The  experiment  was  concluded  before  the  flux  measurements 
attained  a  plateau  indicative  of  an  equilibrium  state  (i.e.,  completion  of  the  drying 
process),  when  an  equilibrium  partition  constant  of  the  dry  soil  could  have  been 
extracted  by  fitting  the  dry  flux  data.  Experiments  performed  at  higher  airflow 
rates  can  provide  more  qualitative  information  for  the  UXO  soil-air  partitioning 
in  the  dry  soils  in  the  time  scales  of  the  experiments  performed. 

A  second  hypothesis  is  that  the  soil  is  dry  enough  to  cause  a  surface  reaction 
to  degrade  the  UXO  compound.  Such  a  loss  via  reaction  can  be  coupled  with 
purely  diffusive  transport  to  model  the  transport  process.  The  gradient  in  the  flux 
data  observed  may  be  representative  of  the  surface  reaction  kinetics  on  dry  soils. 
Since  no  batch  experiments  were  performed  to  evaluate  the  reaction  kinetics  on 
dry  soil,  to  verify  or  evaluate  this  hypothesis  at  this  stage  is  not  possible.  The 
system  is  surmised  to  tend  towards  a  new  equilibrium  state  not  captured  in  the 
time  scale  of  the  experiment.  Previously  reported  oligomerization  reactions 
occur  on  dry  mineral  surfaces  of  soils  with  very  low  organic  carbon.  This 
phenomenon  has  been  observed  by  Karimi-Loftabad,  Pickard,  and  Gray  (1996) 
for  polyaromatic  hydrocarbons  on  dry  soils. 

Figure  8  presents  the  comparison  of  the  experimental  data  for  case  3 
(5-percent  soil  moisture  and  dry  air  over  the  soil)  for  LAAP-C  soil  both  at  20  and 
100  ml/min,  represented  by  symbols  joined  by  continuous  and  broken  lines, 
respectively.  A  significant  observation  between  the  two  sets  of  data  is  that  for 
the  100-ml/min  airflow  rate  case,  the  flux  drops  after  the  first  data  point  and  then 
remains  almost  constant  for  the  duration  of  the  experiment,  while  in  the 
20-ml/min  experiment,  the  flux  drops  continuously,  but  at  a  slower  rate.  In  both 
cases,  the  fluxes  tend  towards  the  same  level  at  the  last  interval  of  experimental 
measurement.  These  observations  suggest  that  the  first  hypothesis,  that  of  a 
transient  stage  of  soil  moisture  resulting  in  a  dynamic  increase  in  the  soil-air 
partition  constant,  is  more  credible  than  the  possibility  of  reaction.  In  conjunc¬ 
tion  with  this,  the  behavior  of  the  experimental  data  in  case  3  at  20  ml/min  is 
more  likely  due  to  slow  continuous  drying  of  the  soil  surface  across  its  length. 
Table  9  lists  the  KSa  fit  values  that  apply  for  the  closest  possible  match  with  the 
experimental  data.  These  values  are  much  higher  than  those  obtained  from  the 
damp  and  wet  soil  cases.  The  KSa  values  listed  for  the  20-ml/min  experiments 
are  listed  along  with  a  k\  value  that  was  used  as  a  fit  parameter  to  fit  the  apparent 
kinetics  in  the  data. 

Case  4, 20-  and  100-ml/min  airflow  experiments.  Figure  9  shows  the 
experimental  data  and  model  fit  for  case  4  with  20-percent  soil  moisture  and 
humid  air  over  the  soil  at  20  and  100  ml/min.  The  transport  of  the  contaminants 
within  the  soil  is  primarily  by  diffusion  in  the  pore  water  followed  by  evapora¬ 
tion  at  the  pore  water-air  interface  at  the  soil  surface.  The  properties  of  the  soil 
are  expected  to  be  constant  over  the  entire  duration  of  the  experiment  since  the 
physical  properties  of  the  soil  do  not  change.  This  is  also  illustrated  by  the 
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Figure  8.  Comparison  of  experimental  fluxes  of  2,4-DNT,  2,6-DNT,  and  1 ,3-DNB  for  case  3  data  with 
LAAP-C  soil  for  20-  and  100-ml/min  airflow  rates 
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Table  9 

Fit  KSa  and  Values  for  Modeling  Experimental  Data  in  Case  3 
from  LAAP-C  Soil 

|  Airflow  Rate 

Ksa ,  L/kg 

fri,  s  1  1 

|  2,4-DNT  || 

20  ml/min 

8.9e+4 

4.7e-6 

100  ml/min 

3.2e+6 

0 

1  2,6-DNT  | 

20  ml/min 

4.8e+4 

3.1e-6 

100  ml/min 

1.3e+6 

0 

1,3-DNB 

20  ml/min 

9.9e+3 

1  .Oe-5 

100  ml/min 

1.3e+6 

0 

experimental  fluxes  that  show  very  little  dynamics.  The  fluxes  at  100  ml/min  are 
higher  than  those  at  20  ml/min  for  almost  all  the  experimental  measurements  in 
this  case.  The  model  fits  were  obtained  by  varying  the  air-water  equilibrium 
constant  values,  KAWi  which  appear  in  the  overall  air-side  mass  transfer  coeffi¬ 
cient.  The  fit  and  the  reported  values  of  KAw  are  listed  in  Table  10.  The  limited 
variation  indicates  the  validity  of  the  model  for  these  conditions. 


Effect  of  soil  temperature  on  flux  from  soils 

Since  temperature  is  an  important  variable  among  the  various  Department  of 
Defense  sites  contaminated  with  explosives  residues,  a  set  of  experiments  was 
conducted  to  study  the  effects  of  soil  temperature  on  flux  to  air.  Figure  10  shows 
the  comparison  of  2,4-DNT  fluxes  from  LAAP-C  soil  samples  at  two  different 
temperatures,  14  °C  and  24  °C  (room  temperature)  for  case  2.  The  flux  was 
0.5 1  ng/cm2/hr  at  14  °C  and  6.63  ng/cm2/hr  at  24  °C  after  24  hr  and  slowly 
decreasing  to  0.1 1  ng/cm2/hr  at  14  °C  and  1.43  ng/cm2/hr  at  24  °C  after  504  hr. 
The  flux  at  14  °C  was  uniformly  lower  than  that  at  24  °C.  Figure  10  also  shows 
the  model  fit  curves  for  the  experimental  data  for  both  temperatures.  The  model 
simulation  for  the  24  °C  case  was  described  earlier  in  this  section.  The  model  fit 
to  the  data  at  14  °C  was  performed  with  the  KsA  obtained  at  24  QC  as  a  starting 
value.  Changing  Ksa  and  minimizing  the  residual  sum  of  square  of  errors 
generated  a  correct  fit.  A  correction  factor  of  ( T2  /Ti)7/4,  where  T  is  temperature 
(°K),  was  used  to  adjust  the  diffusivity  to  the  lower  temperature  (Thibodeaux 
1996).  The  fit  KSA  for  2,4-DNT  on  LAAP-C  soil  at  14  °C  was  107  L/kg  as  com¬ 
pared  to  1.2  x  105  L/kg  at  24  °C.  The  trend  of  higher  partitioning  constant  at 
lower  temperature  was  compared  to  the  data  obtained  for  a  different  compound 
(phenanthrene)  on  a  different  soil,  but  with  comparable  moisture  contents.  The 
heat  of  desorption  from  soil  to  air,  A Hdesorp,  was  calculated  using  the  Clausius- 
Clapeyron  equation 


d(\nP)_  AHdemrp 
dT  RT 2 
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Flux,  ng/cm2/hr  ^ux’  n9/cm2/hr  FIUx,  ng/cm2/hr 


Figure  9.  Comparison  of  experimental  and  model  fit  for  2,4-DNT,  2,6-DNT,  and  1 ,3-DNB  for  case  4  for 
20-  and  100-ml/min  airflow  rates,  LAAP-C  soil,  20-percent  soil  moisture  with  humid  air 
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Table  10 

Fit  Kaw  Values  (or  Henry’s  Law  constant)  to  Model  Experimental 
Data  in  Case  4  from  LAAP-C  Soil 


Parameter 


2,4-DNT 


Reported 


Model  fit,  20  ml/min 


Model  fit,  1 00  ml/min 


2,6-DNT 


Reported 


Model  fit,  20  ml/min 


Model  fit,  100  ml/min 


1,3-DNB 


Reported 


Model  fit,  20  ml/min 


Model  fit,  100  ml/min 


Figure  10.  Effect  of  temperature  on  the  flux  of  2,4-DNT  from  LAAP-C  soil 

(case  2:  20-percent  poisture  content  with  dry  air  passing  over  the  soil) 

where  P  is  the  equilibrium  vapor  pressure  (Pa)  calculated  using  Equation  7 


WSRT 

ksam 
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R  =  universal  gas  constant 
M=  molecular  mass 

A  large  value  for  A Hdesorp  of  1 87  kJ/mol  was  obtained  from  Equation  6  and 
indicated  that  desorption  was  endothermic.  The  heat  of  desorption  for  phenan- 
threne  measured  by  DeSeze  (2000)  was  also  endothermic  (average  value  of 
AHdesorp  for  wet  soil  was  90  ±  4  kJ/mol).  This  suggests  that  the  effects  of  temp¬ 
erature  on  Ksa  and  flux  for  UXO  are  similar  to  those  reported  for  other  types  of 
compounds  and  are  predictable  using  Equation  4  with  appropriate  temperature 
correction  forX^. 


Explosive  Chemical  Signature  Transport 
from  Buried  Sources 

In  these  experiments,  explosive  chemical  signature  transport  from  a  buried 
source  was  evaluated  by  measuring  soil  concentrations  around  and  above  the 
buried  source.  Soil  concentration  measurements  were  acquired  by  sampling 
around  the  buried  source  in  the  experiments  described  earlier.  The  raw  experi¬ 
mental  data  consisted  of  soil  concentration  measurements  at  a  certain  depth  and 
at  a  certain  distance  from  the  central  axis  of  the  container.  The  lateral  sampling 
plane  was  considered  to  be  symmetrical  and,  therefore,  the  measurements  do  not 
indicate  the  radial  direction  of  measurement  about  the  central  axis.  The  raw 
experimental  data  are  plotted  in  Figures  A 1-A8.  These  curves  are  based  on  the 
raw  data  that  were  transformed  by  interpolation  and  smoothing  to  obtain  surface 
plots.  These  plots  are  of  the  soil  concentrations  for  all  four  of  the  explosive 
chemicals  used  in  this  study  (2,4-DNT,  2,6-DNT,  1,4-DNB,  and  1,3,5-TNB)  with 
two  different  soil  conditions  (5  percent  and  20  percent  soil  moisture  content)  and 
two  different  source  depths  (4  cm  and  15.24  cm).  The  Cartesian  coordinates 
X=  0  and  7=  4.0  cm  on  the  figures  corresponds  to  the  source  depth  of  4  cm  in 
the  experimental  data.  The  soil  concentrations  are  higher  nearer  the  source  and 
decrease  away  from  the  source.  In  the  plots  corresponding  to  the  source  depth  at 
15.24  cm,  the  source  coordinates  are  X=  0,  7=  15.24  cm.  The  highest  soil 
concentrations  are  found  closer  to  that  source  at  around  Y-  16  cm  and  decrease 
away  from  the  source.  In  both  cases,  the  soil  concentration  decreases  with 
increasing  absolute  distance  away  from  the  source.  This  is  due  to  the  diffusion 
arising  from  the  concentration  gradients  that  exist  in  the  soil  matrix.  Expectedly, 
the  surface  soil  concentrations  are  higher  in  the  case  where  the  source  is  at  4  cm 
than  when  the  source  is  at  a  depth  of  15.24  cm  from  the  surface.  Also,  the  soil 
concentrations  at  the  surface  immediately  above  the  buried  source  were  higher 
than  on  the  surface  a  few  centimeters  away  from  the  central  axis  of  burial.  Soil 
concentration  measurements  are  listed  and  plotted  to  a  depth  of  about  24  cm  in 
both  cases,  and,  therefore,  no  peaks  are  seen  beyond  7=  20  cm  on  the  plot. 
Deeper  measurements  were  obtained,  but  are  not  shown  here  since  the  interest  is 
primarily  in  the  migration  of  the  chemical  signature  towards  the  surface.  The  5- 
and  20-percent  moisture  data  do  not  differ  qualitatively  or  quantitatively.  This 
was  not  unexpected  because  reaction  rate  tests  for  the  LAAP-C  soil  showed 
similar  results.  The  migration  patterns  of  the  explosive  chemical  seem  similar  in 
both  cases.  The  raw  data  is  also  arranged  in  Tables  A4-A7. 
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To  model  the  transport  of  the  chemical  vapor  signature  away  from  the 
source,  the  Cartesian  coordinates  corresponding  to  the  experimental  measure¬ 
ments  were  converted  to  absolute  distances.  The  quantity,  r,  the  absolute 
distance  from  the  source,  was  calculated  by  the  following  equation: 

r=[(X-X0)2  +  (7-7o)2],/2  (8) 


where 

Xo  and  70  =  source  coordinates 
Xand  7=  sample  coordinates 

Plotting  the  experimental  soil  concentrations  as  a  function  of  r,  one  can  obtain 
the  two-dimensional  plot  to  represent  the  transport  of  the  chemical  from  its 
source.  In  this  analysis,  r  =  0  represents  the  point  at  the  source.  The  chemical 
source  was  assumed  to  be  a  point  source  since  the  plug  containing  the  explosive 
chemicals  was  much  smaller  in  relation  to  the  dimensions  of  the  container. 

FEMLAB,  a  commercially  available  software  from  COMSOL®,  was  initially 
used  to  attempt  to  simulate  a  test  case.  This  program  uses  the  same  diffusion 
model  described  earlier  in  Equation  2  and  follows  a  finite  element  method  to 
solve  for  mobile  phase  concentrations  in  space.  The  program  offers  the  modeling 
capabilities  of  complex  geometries  and  can  generate  surface  plots  that  are 
visually  instructive.  However,  this  program  did  not  provide  a  simple  method  to 
perform  curve  fitting.  An  alternative,  simpler  analytical  model  that  was  more 
amenable  to  fitting  a  parameter  was  also  considered  for  the  specific  task  of 
comparing  and  fitting  the  experimental  data  to  the  model.  This  one-dimensional 
model  tracks  the  concentration  in  soil  as  a  function  of  absolute  distance  away 
from  the  source  and  time.  This  model  uses  the  same  differential  equation  (shown 
in  Equation  2  with  ka  set  to  zero  and  r  replacing  z  in  the  differential  equation),  but 
with  the  following  boundary  conditions:  at  r  =  0  (i.e.,  the  source),  CA=  C°A,  and 
at  r  =  oo  (far  away  from  the  source),  CA  =  0  and  with  the  initial  condition  at  all  r  > 
0,  CA  =  0.  Essentially,  this  assumes  that  at  distances  sufficiently  far  away  from 
the  source,  the  initial  conditions  were  always  maintained.  The  equation  for  the 
soil  concentration  as  a  function  of  time  and  space  derived  from  Equation  2  is 
shown  below  in  Equation  9: 


Ws  (r  ,t)  =  •  CA  ■  erfc 


(9) 


where 

Ws  =  soil  chemical  concentration,  pg/kg 
r  =  absolute  distance  away  from  the  source,  cm 
t  =  time,  hr 

Ksa  =  soil-air  partition  constant,  L/kg 
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C°A  =  soil  vapor  concentration  at  the  source,  pg/cm3 

a  =DeJ/Rf,  where  Dejj  and  7?/ are  the  effective  diffusivity  (cm2/hr)  and  the 
retardation  factor,  respectively 

For  the  5-percent  soil  moisture  case,  the  air-filled  porosity  of  the  soil  was 
assumed  to  be  constant  since  there  could  not  be  any  significant  evaporation.  The 
loss  of  water  vapor  to  the  headspace  in  the  containers  was  considered  negligible 
and  not  to  affect  soil  porosity  significantly.  For  the  5-percent  soil,  the  soil  air 
porosity  was  set  at  0.23,  and  vapor-phase  diffusion  was  considered  to  be  the 
primary  mechanism  of  chemical  transport.  The  buried  source  was  assumed  to  be 
a  constant  source  of  saturated  vapor  of  each  of  the  chemicals  under  considera¬ 
tion.  The  vapor-phase  concentration  at  source  was  calculated  from  the  reported 
vapor  pressure  and  at  room  temperature.  These  computed  values  were  0.0021, 
0.0055,  and  0.0012  pg/cm3  for  2,4-DNT,  2,6-DNT,  and  1,3-DNB,  respectively. 

In  the  20-percent  soil  moisture  case,  diffusion  in  pore  water  was  considered 
the  primary  transport  mechanism  since  the  soil  pore  spaces  were  filled  with 
water.  The  water- filled  porosity  was  set  at  0.31.  In  Equation  7,  Ksa  was  replaced 
by  Ksw,  C°A  was  the  solubility  of  the  chemical  in  water,  DCfj  was  computed  using 

the  diffusivity  of  the  chemical  in  water,  and  Rf  was  calculated  using  Ksw •  The 
saturated  aqueous  concentrations  were  280,  208,  and  533  pg/cm3  for  2,4-DNT, 
2,6-DNT,  and  1,3-DNB,  respectively.  The  same  model  was  applied  to  the  data 
from  two  different  burial  depths.  For  the  5-percent  soil  moisture  data,  a  best-fit 
model  was  applied  to  the  experimental  data  by  varying  Ksa •  Figure  1 1  shows  the 
comparison  between  the  model  fit  and  the  experimental  data.  The  model  follows 
the  general  trend  shown  by  the  experimental  soil  concentrations,  which  decrease 
away  from  the  source.  Table  1 1  lists  the  Ksa  for  the  model  fit  curves  shown  in 
Figure  1 1  and  also  compares  these  values  to  those  obtained  earlier.  The  Ksa 
values  vary  over  a  range  of  two  orders  of  magnitude.  In  previous  experiments, 
Jenkins  et  al.  (2000)  measured  KSa  values  for  explosive  chemicals  and  found  that 
values  varied  from  3.9e+4  to  1  .Oe+6  L/kg  for  2,4-DNT,  2.6e+4  to  4.1e+5  L/kg 
for  2,6-DNT,  and  5.34e+4  to  2.8e+6  L/kg  for  1,3-DNB.  The  values  measured  by 
Jenkins  et  al.  (2000)  were  for  a  wet  soil  with  moisture  contents  greater  than 
10  percent  and  are  comparable  to  the  range  of  fit  KSa  values  obtained  from  the 
current  study.  Figure  12  shows  the  comparison  between  the  model  and  experi¬ 
mental  soil  concentrations  for  the  20-percent  soil  moisture  case.  The  model 
curves  were  obtained  by  using  the  reported  aqueous  solubility  and  the  measured 
Ksw  values.  The  experimental  and  the  model  curves  match  very  well.  The  model 
is  a  good  predictor  of  the  transport  of  the  UXO  chemicals  from  a  buried 
munitions  source  when  the  source  concentration  is  known. 

Table  12  lists  the  comparisons  of  the  air  concentrations  calculated  from  the 
experimental  and  model  fit  soil  concentrations.  For  the  5-percent  soil  moisture 
cases,  the  experimental  air  concentrations  were  obtained  by  dividing  the  experi¬ 
mental  soil  concentrations  by  the  fit  Ksa  values  shown  in  Table  1 1 .  For  the 
20-percent  soil  moisture  cases,  the  estimated  KSa  from  the  reported  Ksw  and  KAW 
was  used  (also  shown  in  Table  1 1).  For  the  5-percent  soil  moisture  case,  the 
model  directly  provides  the  vapor-phase  concentrations.  For  the  20-percent  soil 
moisture  case,  the  model  solves  for  aqueous  concentrations;  therefore,  at  the 
surface  the  aqueous  concentrations  were  multiplied  by  KAw  to  provide  the  air 
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Concentration  in  Soil,  pg/pg  Concentration  in  Soil,  pg/pg  Concentration  in  Soil,  pg/pg 


Figure  11.  Comparison  of  experimental  and  modeled  soil  concentrations  of  2,4-DNT,  2,6-DNT,  and 
1 ,3-DNB  from  a  buried  source,  soil  moisture  of  5  percent  with  vapor-phase  diffusion 
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Table  11 

Adjustable  Parameters,  KSa  (L/kg),  for  Comparison  with  Laboratory 
Flux  Experiments  and  Buried  Chemical  Experiments  with  LAAP-C 
Soil 

||  Parameter 

Case  1  Case  2 

I  2,4-DNT 

Estimated 

8.9e+4 

8.9e+4 

Model  fit  -  20  ml/min 

1.2e+5 

9.7e+4 

Model  fit  - 100  ml/min 

1.6e+5 

2.9e+4 

Model  fit  -  Buried  (4  cm)  (5%  soil  moisture) 

1.1e+3 

- 

Model  fit  -  Buried  (15  cm)  (5%  soil  moisture) 

1.4e+3 

- 

|  2,6-DNT 

Estimated 

1 .3e+5 

1 .3e+5 

Model  fit  -  20  ml/min 

6.4e+3 

2.5e+4 

Model  fit  - 100  ml/min 

1.5e+4 

3.2e+4 

Model  fit  -  Buried  (4  cm)  (5%  soil  moisture) 

2.1e+3 

- 

Model  fit  -  Buried  (15  cm)  (5%  soil  moisture) 

3.0e+3 

- 

|  1,3-DNB 

Estimated 

9.9e+3 

9.9e+3 

Model  fit  -  20  ml/min 

6.1e+4 

6.9e+4 

Model  fit  - 100  ml/min 

1 .1e+5 

2.0e+4 

Model  fit  -  Buried  (4  cm)  (5%  soil  moisture) 

2.0e+3 

- 

Model  fit  -  Buried  (15  cm)  (5%  soil  moisture) 

3.4e+3 

- 

concentrations.  The  predicted  air  concentrations  are  of  the  same  order  of  magni¬ 
tude,  or  slightly  higher  than,  all  of  the  measured  air  concentrations,  except  for  the 
5-percent  soil  moisture  case  with  the  source  at  4  cm  from  the  surface.  The  mea¬ 
sured  concentrations  were  higher  than  the  predicted  concentration.  In  general, 
the  estimation  of  KSa  seemed  to  be  the  most  significant  factor  in  predicting  the 
soil  and  air  concentrations  around  a  buried  explosive  source.  The  estimation  of 
Ksa  from  Ksw  and  Kwa  may  not  be  very  appropriate.  The  Ksa  must  be  measured 
empirically  for  the  desired  soil  by  varying  parameters  such  as  moisture  content  or 
temperature.  A  gas  saturation  method  has  been  used  successfully  in  the  past  to 
measure  the  soil-air  equilibrium  partition  constants  for  several  polyaromatic 
hydrocarbons  (DeSeze  2000).  In  this  method,  clean  air  with  the  desired  relative 
humidity  is  passed  through  a  soil  column  of  sufficient  length  to  facilitate  equilib¬ 
rium.  The  clean  air  exiting  the  column  is  sent  through  a  plug  of  adsorbent 
material.  Both  the  air  and  soil  concentrations  are  analyzed  to  obtain  Ksa-  An 
alternative  method,  described  by  Jenkins  et  al.  (2000),  measured  Ksa  by  placing  a 
plug  of  adsorbent  material  over  contaminated  soil  and  analyzing  the  plug  and  the 
soil  after  equilibrium  is  achieved. 


Application  of  the  Model  to  Soil  Samples 
Collected  Near  Buried  Landmines 

The  study  conducted  by  Jenkins  et  al.  (2000)  to  measure  explosive-related 
signatures  in  soil  samples  near  buried  landmines  provided  a  data  set  that  could  be 
used  to  verify  the  model.  This  study  was  conducted  at  Fort  Leonard  Wood,  MO. 
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a)  2,4-DNT  -  Source  4  cm  deep 


Experimenfa! 
•  Model  Fit 


10  20 
Distance  from  Source,  cm 


b)  2,4-DNT  -  Source  15  cm  deep 


10  20 
Distance  from  Source,  cm 


c)  2,6-DNT  -  Source  4  cm  deep 


Experimental 
■  Model  Fit 
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10  20 
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d)  2,6-DNT  -  Source  15  cm  deep 


10  20 
Distance  from  Source,  cm 


e)  1 ,3-DNB  -  Source  4  cm  deep 


f)  1 ,3-DNB  -  Source  1 5  cm  deep 


10  20 
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10  20 
Distance  from  Source,  cm 


Figure  12.  Comparison  of  experimental  and  modeled  soil  concentrations  of  2,4-DNT,  2,6-DNT,  and 
1,3-DNB  from  a  buried  source,  soil  moisture  of  20  percent  with  aqueous-phase  diffusion 
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Table  12 

Comparison  of  Air  Concentrations  Calculated  from  the  Buried 
Source  Experimental  Data  and  Corresponding  Model 

Simulations 

Soil  Moisture,  % 

Source  Depth,  cm 

Air  Concentration,  pg/cm3 

Experimental  Model 

2,4-DNT 

5 

4 

3.17e-2 

1.33e-3 

5 

15 

3.59e-5 

8.12e-5 

20 

4 

7.97e-5 

8.49e-4 

20 

15 

6.34e-7 

3.21e-6 

2,6-DNT  1 

5 

4 

6.38e-3 

2.8e-3 

5 

15 

1.91e-4 

1.34e-5 

20 

4 

2.77e-5 

5.25e-4 

20 

15 

1 .34e-6 

4.06e-6 

1,3-DNB  I) 

5 

4 

4.01  e-2 

6.84e-4 

5 

15 

9.8e-5 

1.81e-5 

20 

4 

1.14e-3 

9.01e-3 

20 

15 

2.33e-5 

2.83e-4 

Several  types  of  landmines  were  buried  at  different  depths  in  soil,  and  core  sam¬ 
pling  was  performed  at  several  time  intervals  up  to  472  days  after  the  burial.  The 
soil  concentrations  were  measured  in  the  collected  samples  and  were  tabulated  as 
a  function  of  the  mine  type,  collection  date,  and  collection  location.  Samples 
were  collected  all  around  the  mine  and  also  below  it  after  the  mine  was  removed. 
One  of  the  observations  made  in  the  Jenkins  study  was  that  high  concentrations 
of  explosive  chemicals  were  found  under  the  mine,  possibly  due  to  leaching. 

Researchers  chose  to  compare  the  model  simulations  with  two  of  the  data 
sets  from  one  mine  type  (TMA-5)  at  two  different  dates  of  sampling  (August 
1998  and  April  1999)  that  were  41  and  270  days  after  the  mines  were  buried. 

The  attempt  was  to  compare  the  model  simulations  with  the  measured  concen¬ 
trations  above  the  depth  of  the  mine  (the  TMA-5  mines  were  buried  at  a  depth  of 
8.5  cm  below  the  soil  surface).  Since  leaching  was  reported  below  the  depth  of 
the  mine,  vapor  transport  likely  occurred  predominantly  in  an  upward  direction 
from  the  mine.  The  simple  model  described  in  the  earlier  section  was  used  to 
simulate  soil  concentrations  at  a  distance  away  from  the  source.  The  mine  was 
not  a  symmetric  source,  but  was  31  cm  long,  14  cm  wide,  and  1 1  cm  high.  The 
distance  from  the  source  was  computed  from  the  nearest  edge  of  the  mine.  Infor¬ 
mation  about  the  lateral  positions  of  the  sampling  cores  and  sampling  depths  in 
Jenkins  et  al.  (2000)  was  used  to  obtain  the  coordinates  of  the  samples.  The 
mean  depth  of  a  sample  in  a  core  was  assumed  to  be  the  depth  of  the  sample  (i.e., 
if  a  sample  represented  5  to  10  cm  of  a  core,  the  depth,  or  Y-coordinate,  of  the 
sample  was  set  at  7.5  cm).  The  porosity  of  the  soil  was  computed  using  the  par¬ 
ticle  size  distribution  information  and  was  set  at  0.47.  The  bulk  density  informa¬ 
tion  was  not  found  in  Jenkins  et  al.  (2000)  and  hence  a  value  of  1 .43  g/cm3 
(based  on  LAAP-C  soil)  was  assumed.  The  model  also  assumed  source  vapor- 
phase  concentrations  based  on  the  reported  vapor  pressures  at  room  temperatures. 
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The  initial  KSa  values  used  in  the  simulation  were  the  measured  values  reported 
in  their  study.  The  simulations  using  these  parameters  overpredicted  the  soil 
concentrations  by  several  orders  of  magnitude  within  the  first  2  to  5  cm  from  the 
source  and  underpredicted  the  soil  concentrations  for  distances  beyond  5  cm.  In 
other  words,  the  orientation  and  curvature  of  the  model  curves  cut  across  the 
experimental  data  at  almost  right  angles.  Adjusting  the  Ksa  values  alone  changed 
the  curvature  of  the  model  simulation  curves,  but  the  magnitude  was  a  few  orders 
higher  than  the  experimental  data.  Adjusting  the  vapor-phase  concentration  at 
the  source  was  also  necessary  to  obtain  good  fit  between  the  model  and  the  data. 
Figures  13  and  14  show  the  comparison  between  the  model  fit  and  the  experi¬ 
mental  data.  The  adjusted  values  for  2,4-DNT,  2,6-DNT,  1,3-DNB,  and  TNT  are 
listed  in  Table  13  along  with  the  reported  values.  The  variation  in  the  Ksa  does 
not  seem  out  of  place  considering  the  variation  observed  in  the  measured  values 
and  the  model  fit  in  the  laboratory  flux  experiments  reported  in  previous  sections. 
However,  the  saturated  vapor  concentrations  input  in  the  model  fit  are  a  few 
orders  of  magnitude  different  from  reported  values.  The  median  air  concentra¬ 
tion  of  2,4-DNT  at  the  soil-air  interface  on  top  of  a  TMA-5  mine  buried  8.5  cm 
was  estimated  on  the  basis  of  measured  Ksa  to  be  1.5  x  10'7  pg/cm3.  As  a  com¬ 
parison,  the  surface  concentrations  estimated  from  the  laboratory  experiment 
with  a  source  buried  4  cm  and  15  cm  under  the  surface  were  0.03 17  and 
7.9  x  10‘5  pg/cm3,  respectively. 

A  significant  difference  between  the  studies  of  buried  sources  conducted  in 
the  laboratory  and  the  field  study  data  used  to  test  the  model  is  the  structure  and 
strength  of  the  explosive  chemical  source.  In  the  laboratory  study,  the  source 
was  a  small  plug  that  was  on  the  order  of  a  few  centimeters  and  was  in  direct 
contact  with  the  soil.  Therefore,  the  transport  processes  in  the  porous  media 
could  be  easily  visualized  and  modeled.  In  the  case  of  munitions,  chemical 
vapors  are  generated  inside  the  casing  and  travel  to  the  surrounding  soil  through 
the  casing  unless  the  casing  is  breached.  This  presents  an  additional  obstacle  for 
vapor-phase  transport  to  soil  that  is  not  easily  overcome.  If  a  casing  is  breached, 
it  cannot  be  considered  a  symmetric  source  and  may  result  in  asymmetrical 
movement  of  the  signature  in  soil.  The  largest  obstacle,  however,  to  both  model¬ 
ing  and  detection  of  chemical  signatures  may  be  the  lack  of  a  consistent  and 
known  signature  from  UXO.  This  is  borne  out  in  the  study  of  Phelan,  Rodacy, 
and  Barnett  (2001)  who  found  that  the  chemical  signature  emitted  from  simulated 
unexploded  ordnance  was  insufficient  for  use  in  discrimination  of  live  explosive- 
containing  items.  A  unique,  persistent,  distinguishing  trace  chemical  signature 
emitted  from  UXO  was  not  found  under  all  test  and  field  conditions.  A 
necessary  precursor  to  a  good  predictive  model  is  a  known  source  term.  The 
model  developed  in  this  study  can  predict  with  good  accuracy  the  transport  of 
explosives  signatures  and  resulting  soil  concentrations  at  a  distance  from  the 
source  when  the  source  is  in  contact  with  the  soil.  This  will  be  useful  in  both 
surface  and  vadose  zone  soils. 
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Figure  13.  Comparison  of  experimental  (data  from  Jenkins  et  al.  (2000))  and  model  fit  soil  concentrations 
of  2,4-DNT,  2,6-DNT,  1,3-DNB,  and  TNT  from  buried  TMA-5  type  landmine,  Fort  Leonard 
Wood  site,  41  days  following  burial 
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Soil  Concentration,  jxg/kg  Soil  Concentration,  pg/kg 


Figure  14.  Comparison  of  experimental  (data  from  Jenkins  et  al.  (2000))  and  model  fit  soil  concentrations 
of  2,4-DNT,  2,6-DNT,  1,3-DNB,  and  TNT  from  buried  TMA-5  type  landmine,  Fort  Leonard 
Wood  site,  270  days  following  burial 
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Table  13 

Model  Fit  Parameters  for  the  Comparison  of  Landmine  Data  with 
Model  Simulations 

Parameter 

2,4-DNT 

2,6-DNT 

1,3-DNB 

TNT 

K$a ,  L/kg,  estimated 

8.6e+5 

1.8e+5 

4.6e+5 

9.6e+6 

Ksa ,  L/kg,  fit 

5.1e+3 

1.2e+3 

4.1e+3 

5.3e+3 

Vapor  concentration, 
pg/cm5,  estimated 

2.1  e-3 

5.5e-3 

1 ,2e-3 

6.7e-5 

Vapor  concentration, 

^ig/cm3,  fit 

3.2e-5 

5.1e-6 

5.1e-6 

5.2e-6 
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Increasing  temperature  and  moisture  in  a  surface  soil  decreased  the  stability 
of  1,3,5-TNB,  1,3-DNB,  TNT,  2,4-DNT,  and  2,6-DNT.  In  contrast,  concentra¬ 
tions  of  explosives  in  the  aquifer  soils  remained  unchanged  regardless  of  temp¬ 
erature  and  moisture.  Explosives  signatures  in  surface  soils  will  therefore 
undergo  greater  degradation  than  in  subsurface  soils,  making  surface  detection  of 
signatures  more  difficult. 

The  moisture  content  and  temperature  of  soils  affect  the  flux  of  explosives 
chemical  signatures  from  soil  surfaces.  The  drier  a  soil  becomes,  the  more  the 
explosives  signature  flux  decreases,  probably  due  to  increased  partitioning  on  the 
soil.  Fluxes  of  explosives  signatures  are  relatively  constant  for  soils  with  at  least 
5-percent  moisture  and  that  do  not  dry  out.  Increased  temperature  results  in  order 
of  magnitude  increases  in  explosives  signature  fluxes  from  soil.  Explosive  signa¬ 
tures  from  soils  can  be  expected  to  increase  as  the  soils  become  wetter  and 
temperatures  rise. 

The  compounds  2,6-DNT,  1,3-DNB,  and  2,4-DNT  were  most  likely  to  pro¬ 
vide  chemical  signatures  from  soils.  Of  these  compounds,  2,6-DNT  generally 
provided  the  strongest  signature.  However,  even  the  strongest  2,6-DNT  signa¬ 
tures  were  relatively  low,  ranging  from  0.1  to  <  100  ng/cm2/hr,  depending  upon 
the  soil.  These  low  signatures  mean  that  sensors  will  need  to  be  very  sensitive, 
which  may  result  in  many  false  alarms  because  of  explosives  residues  not 
associated  with  UXO  on  firing  ranges.  Low  fluxes  of  even  the  most  volatile 
compounds  from  explosives  indicate  that  this  environmental  loss  pathway  is 
minimal  and  should  not  result  in  significant  environmental  exposures. 

The  model  and  process  descriptors  developed  for  prediction  of  explosives 
signature  movement  in  soils  accurately  predicted  signature  movement  in  a  soil 
with  20-percent  moisture  content  and  explosives  sources  buried  at  various 
depths.  For  drier  soils,  site-specific  measurements  of  the  soil-air  partitioning 
coefficient  may  be  needed  pending  development  of  improved  Henry’s  Law 
constants.  Verification  of  the  model  using  an  independent  data  set  for  explosives 
signatures  around  landmines  showed  that  adjusting  the  vapor-phase  concentration 
at  the  source  was  necessary  to  obtain  a  good  model  fit.  This  was  needed  because 
the  casing  of  munitions  presents  a  barrier  to  signature  movement  that  is  not  easily 
overcome  and  results  in  an  inconsistent  and  variable  vapor  emission.  The  model 
can  accurately  predict  explosives  signature  movement  to  the  surface  where 
chemical  detection  can  occur  when  the  source  strength  is  known.  The  model  can 
also  predict  explosives  signature  movement  and  corresponding  accumulation  of 
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explosives  concentrations  in  vadose  zone  soils.  The  model  verification  phase  of 
this  study  also  revealed  the  need  for  improved  Henry’s  Law  constants  that  reflect 
field  conditions. 
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Figure  A1 .  Experimental  soil  concentrations  of  2,4-DNT  as  a  function  of  position  with  respect  to  source  at  4-cm  depth 
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Figure  A2.  Experimental  soil  concentrations  of  2,4-DNT  as  a  function  of  position  with  respect  to  source  at  15-cm  depth 


Source  Coordinates:  X  =  0  Y  =  4.064 
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Figure  A3.  Experimental  soil  concentrations  of  2,6-DNT  as  a  function  of  position  with  respect  to  source  at  4-cm  depth 
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Figure  A4.  Experimental  soil  concentrations  of  2,6-DNT  as  a  function  of  position  with  respect  to  source  at  1 5-cm  depth 
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Figure  A5.  Experimental  soil  concentrations  of  1 ,3-DNB  as  a  function  of  position  with  respect  to  source  at  4-cm  depth 
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Figure  A6.  Experimental  soil  concentrations  of  1 ,3-DNB  as  a  function  of  position  with  respect  to  source  at  15-cm  depth 
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Figure  A7.  Experimental  soil  concentrations  of  1 ,3,5-TNB  as  a  function  of  position  with  respect  to  source  at  4-cm  depth 


Source  Coordinates 
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Figure  A8.  Experimental  soil  concentrations  of  1,3,5-TNB  as  a  function  of  position  with  respect  to  source  at  15-cm  depth 


Table  A1 
Initial  an< 
Temperat 

1 

i  Final  Concentration  (mg/kg)  in  Soils  as  a  Function  of 
ture  and  Moisture 

Compound 

Initial 

Concentration 

14  °C 

24  °C  | 

5%  moisture 

20%  moisture 

5%  moisture 

20%  moisture  [ 

LAAP-C  soil  following  28  days  of  incubation  | 

1 ,3,5-TNB 

7.56  ±  0.02 

7.75  +  0.16 

7.63  +  0.09 

7.12  +  0.16 

7.29  +  0.21 

1,3-DNB 

9.98  +  0.12 

10.85  +  0.35 

10.45  +  0.25 

9.84  +  0.16 

9.63  +  0.37 

TNT 

9.97  +  0.13 

10.55  +  0.35 

10.2  +  0.1 

9.65  +  0.18 

9.60  +  0.33 

2,6-DNT 

10.4  +  0.1 

10.23  +  0.57 

10.03  +  0.27 

9.76  +  0.11 

8.99  +  0.48 

2,4-DNT 

10.35  +  0.15 

9.97  +  0.54 

9.84  +  0.04 

9.36  +  0.16 

9.55  +  0.22 

I  LAAP-D  soil  following  28  days  of  incubation  | 

1,3,5-TNB 

8.17  +  0.19 

7.87  +  0.26 

8.64  +  0.06 

7.79  +  0.44 

8.48  +  0.44 

1,3-DNB 

10.4  +  0.1 

10.7  +  0.3 

11.85  +  0.05 

10.75  +  0.35 

11.6  +  0.6 

TNT 

10.8  +  0.3 

10.6  +  0.4 

11.65  +  0.05 

10.73  +  0.37 

11.4  +  0.6 

2,6-DNT 

11  +0.3 

11+0.1 

11.6  +  0.2 

10.7  +  0.5 

11.85  +  0.75  j 

2,4-DNT 

11  +0.2 

10.55  +  0.15 

11.5  +  0.1 

10.37  +  0.54 

11.6  +  0.7  | 

j  Yokena  cay  soil  following  14  days  of  incubation  t\ 

1,3,5-TNB 

7.51+0.15 

6.48  +  0.07 

<0.1 

6.09  +  0.11 

<0.1 

1,3-DNB 

9.76  +  0.04 

9.70  +  0.02 

<0.1 

8.99  +  0.04 

<0.1 

TNT 

9.65  +  0.13 

8.45  +  0.05 

0.46  +  0.005 

7.57  +  0.13 

0.20  +  0 

2,6-DNT 

10.15  +  0.05 

9.43  +  0.01 

7.96  +  0 

8.89  +  0.01 

0.26  +  0.02 

2,4-DNT 

10.0  +  0 

10.1  +0 

2.07  +  0.005 

9.65  +  0.5 

0.24  +  0.02 
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Table  A2 

Mass  (pg)  of  Explosives  from  Buried  Explosives  at  5-Percent  Soil 
Moisture  and  4-cm  Burial  Depth 


Horizontal/Vertical 

Distance  from  Source,  cm 

1,3,5-TNB 

TNT 

2,4-DNT 

0  cm  / 10.2  cm  deep 

227 

153 

168 

183 

131 

0  cm  / 15.3  cm  deep 

214 

56.8 

54.8 

88.5 

67.8 

0  cm  /  20.4  cm  deep 

18.7 

26.3 

10.8 

25.2 

15.4 

0  cm  /  23.8  cm  deep 

1.53 

1.56 

1.20 

2.65 

1.33 

2.5  cm  /  0  cm  (top  cm) 

11.7 

488 

4.32 

108 

282 

2.5  cm  /  3.4  cm  deep 

233 

163 

40.6 

39.9 

112 

2.5  cm  /  8.5  cm  deep 

280 

229 

135.7 

192 

153 

6.5  cm  / 1.7  cm  deep 

0.142 

40.7 

0.063 

20.9 

19.5 

6.5  cm /11.9cm  deep 

0.071 

1.84 

0.047 

3.30 

0.736 

6.5  cm  / 17  cm  deep 

0.066 

1.73 

0.038 

1.97 

0.694 

6.5  cm  /  25.5  cm  deep 

0.062 

1.43 

0.050 

0.791 

9.5  cm  /  0  cm  (top  cm) 

0.084 

5.97 

0.025 

11.4 

1.34 

9.5  cm  /  6.8  cm  deep 

0.062 

0.414 

0.023 

3.29 

0.078 

9.5  cm /23.8  cm  deep 

0.079 

0.517 

0.045 

1.48 

0.094 

12  cm  / 1.7  cm  deep 

0.074 

0.588 

0.022 

2.75 

0.137 

12  cm  / 13.6  cm  deep 

0.078 

0.110 

0.034 

0.525 

0.046 

16  cm  /  5.1  cm  deep 

0.062 

0.149 

0.023 

1.12 

0.062 

16  cm  /  8.5  cm  deep 

0.077 

0.155 

0.034 

0.677 

0.062 

16  cm/ 20.4  cm  deep 

0.064 

0.061 

0.021 

0.165 

0.043 

16  cm /23.8  cm  deep 

0.057 

0.041 

0.044 

0.184 

0.026 
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Table  A3 

Mass  (pg)  of  Explosh 
Moisture  and  15-cm  1 

ves  from  B 
Burial  Dep 

uried  Explosives  at  5-Percent  Soil 
th 

Horizontal/Vertical 

Distance  from  Source,  cm 

1,3,5-TNB 

1,3-DNB 

TNT 

0  cm  /  0  cm  (top  cm) 

13.9 

1.12 

0.290 

3.26 

0.287 

0  cm  /  20.4  cm  deep 

214 

123 

158 

212 

153 

2.5  cm  /  5.1  cm  deep 

1.21 

4.52 

0.308 

10.8 

0.942 

2.5  cm  /  23.8  cm  deep 

13.0 

10.4 

3.59 

12.2 

11.3 

5.5  cm  / 1.7  cm  deep 

0.103 

0.302 

0.089 

2.78 

0.122 

5.5  cm  / 11.9  cm  deep 

0.112 

28.8 

0.108 

29.3 

18.6 

5.5  cm /17  cm  deep 

0.088 

13.4 

0.080 

23.1 

9.11 

5.5  cm  /  27.2  cm  deep 

0.076 

4.49 

0.044 

8.04 

3.80 

8.5  cm  /  6.8  cm  deep 

0.063 

0.521 

0.040 

6.99 

0.094 

8.5  cm  /  23.8  cm  deep 

0.056 

0.816 

0.028 

3.99 

0.123 

11.5  cm  /  0  cm  (top  cm) 

0.111 

0.144 

0.022 

1.15 

0.153 

11.5  cm  / 13.6  cm  deep 

0.056 

0.124 

0.024 

2.37 

0.057 

15.5  cm  /  5.5  cm  deep 

0.055 

0.079 

0.022 

0.714 

0.059 

15.5  cm  / 15.3  cm  deep 

0.054 

0.030 

0.031 

0.970 

0.046 

15.5  cm  /  20.4  cm  deep 

0.058 

0.112 

0.027 

1.21 

0.076 

15.5  cm  /  27.2  cm  deep 

0.056 

0.017 

0.019 

0.314 

0.044 
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Table  A4 

Mass  (pg)  of  Explosives  from  Buried  Explosives  at  20-Percent  Soil 
Moisture  and  4-cm  Burial  Depth 


Horizontal/Vertical 

Distance  from  Source,  cm 

RBI 

1,3-DNB 

TNT 

IHHHURI 

0  cm  / 10.2  cm  deep 

45.6 

99.5 

69.8 

46.1 

67.0 

0  cm  / 15.3  cm  deep 

10.4 

34.3 

36.0 

32.4 

37.9 

0  cm  /  20.4  cm  deep 

3.31 

15.6 

18.4 

17.6 

21.3 

0  cm  /  27.2  cm  deep 

6.35 

45.8 

26.8 

36.3 

44.1 

2.5  cm  /  0  cm  (top  cm) 

42.2 

53.3 

46.4 

16.7 

33.5 

2.5  cm  / 1.7  cm  deep 

79.0 

84.5 

76.5 

31.2 

52.4 

2.5  cm  /  3.4  cm  deep 

111 

101 

94.1 

48.2 

77.5 

2.5  cm  /  5.1  cm  deep 

162 

147 

333 

207 

185 

2.5  cm  /  8.5  cm  deep 

130 

126 

80.9 

39.5 

79.7 

4.5  cm  / 1.7  cm  deep 

14.3 

39.5 

12.9 

12.8 

20.4 

4.5  cm  / 11.9  cm  deep 

0.141 

34.1 

4.01 

8.74 

9.32 

4.5  cm  / 17  cm  deep 

<0.047 

11.6 

1.92 

4.63 

4.56 

7  cm  /  0  cm  (top  cm) 

0.087 

40.6 

5.41 

9.08 

11.4 

7  cm  /  6.8  cm  deep 

11.5 

185 

35.3 

49.6 

82.4 

7  cm  /  23.8  cm  deep 

<0.047 

0.059 

<0.047 

0.149 

<0.047 

9.5  cm  / 1 .7  cm  deep 

<0.047 

10.9 

0.284 

3.24 

2.51 

9.5  cm  / 13.6  cm  deep 

<0.047 

2.21 

<0.047 

1.21 

0.158 

12  cm  / 10.2  cm  deep 

<0.047 

1.38 

<0.047 

1.04 

<0.047 

12  cm  /  20.4  cm  deep 

<0.047 

0.211 

<0.047 

0.248 

<0.047 

14.5  cm  /  0  cm  (top  cm) 

<0.047 

2.56 

<0.047 

2.02 

0.789 

14.5  cm  /  5.1  cm  deep 

<0.047 

0.766 

<0.047 

0.911 

<0.047 

14.5  cm  /  28.9  cm  deep 

<0.047 

<0.047 

<0.047 

<0.047 

<0.047 

17  cm  /  8.5  cm  deep 

<0.047 

0.217 

<0.047 

0.583 

<0.047 

17  cm  / 17  cm  deep 

<0.047 

0.073 

<0.047 

0.107 

<0.047 
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Table  A5 

Mass  (pg)  of  Explosives  from  Buried  Explosives  at  20-Percent  Soil 
Moisture  and  15-cm  Burial  Depth 


Horizontal/Vertical 

Distance  from  Source,  cm 

1,3,5-TNB 

TNT 

0  cm  /  0  cm  (top  cm) 

<0.047 

0.716 

<0.047 

0.533 

0.175 

0  cm  /  8.5  cm  deep 

12.8 

29.3 

6.90 

13.1 

18.9 

0  cm  /  20.4  cm  deep 

142 

70.3 

121 

82.6 

105 

0  cm  /  25.5  cm  deep 

114 

49.0 

79.8 

77.6 

76.3 

2.5  cm  /  5.1  cm  deep 

3.14 

6.83 

0.924 

3.67 

3.53 

2.5  cm /23.8  cm  deep 

0.158 

5.07 

0.223 

2.80 

1.42 

5.5  cm  / 1.7  cm  deep 

<0.047 

1.09 

0.039 

1.10 

0.193 

5.5  cm  / 11.9  cm  deep 

15.2 

55.5 

16.1 

27.9 

45.6 

5.5  cm  / 17  cm  deep 

0.462 

20.7 

7.44 

16.2 

23.5 

5.5  cm  /  27.2  cm  deep 

<0.047 

0.252 

0.078 

0.446 

0.063 

7.5  cm  /  6.8  cm  deep 

<0.047 

7.56 

0.140 

4.28 

2.34 

7.5  cm  /  23.8  cm  deep 

<0.047 

1.19 

0.052 

1.05 

0.093 

1 0  cm  /  0  cm  (top  cm) 

<0.047 

0.151 

<0.047 

0.346 

<0.047 

10  cm  / 13.6  cm  deep 

<0.047 

6.39 

0.778 

3.78 

1.69 

13  cm  /  20.4  cm  deep 

0.104 

0.789 

0.163 

0.701 

0.223 

15  cm  /  3.4  cm  deep 

0.049 

0.039 

0.019 

0.229 

0.018 

15  cm  / 10.2  cm  deep 

0.061 

0.201 

0.049 

0.451 

0.030 

15  cm  /  28.9  cm  deep 

0.073 

<0.047 

0.111 

0.049 

<0.047 

17.5  cm  / 15.3  cm  deep 

<0.047 

0.111 

<0.047 

0.469 

0.021 

17.5  cm  /23.8  cm  deep 

<0.047 

<0.047 

<0.047 

0.056 

<0.047 
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Appendix  B 
Notation 


Ac  Sbil-air  interface  area  of  the  chamber 

C  Chemical  concentration  of  reacting  substance 

CA  Concentration  of  chemical  in  pore  air 

Co  Concentration  of  the  reactant  at  time  0 

C°A  Initial  equilibrium  concentration  of  chemical  in  pore  air  (kg/m3) 

C°r  Equilibrium  vapor-phase  concentration  of  the  chemical  in  soil 

d  Length  of  the  airflow  path 

Deff  Effective  diffusivity  (m2/s) 

Da  Diffusivity  of  contaminant  A  in  air  (m2/s) 

H(  Henry’s  Law  constant 

ka  Air-side  mass  transfer  coefficient  (m/s) 

kaCh  Advective  velocity  (cm3/cm2/hr) 

k\  First-order  constant  for  disappearance  rate  of  chemical  in  soil  (1/s) 

Kaw  Air-water  partition  constant  or  Henry’s  Law  constant 

(dimensionless) 

Ksa  Soil-air  partition  constant  (m3/kg) 

Ksw  Soil-water  partition  constant  (L/kg) 

M  Molecular  mass 

Na  Flux  rates 

P  Equilibrium  vapor  pressure,  Pa 

r  Absolute  distance  away  from  a  chemical  source  (cm) 

R  Universal  gas  constant 

Re  Reynolds  number 

Rf  Retardation  factor  (dimensionless) 

Sc  Schmidt  number 


B1 
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T 

v 

Ws(r,f) 

W° 


*>y 


X0,Y0 

z 

AHdesorp 

Am 

At 

la 

Y* 

It 

Pb 

V 


Temperature  (°K) 

Average  flow  velocity 

Soil  concentration  or  loading  of  contaminant  (kg/kg) 

Initial  soil  concentration  or  loading  of  contaminant  (kg/kg) 
Sample  coordinates 
Source  coordinates 
Vertical  distance  (cm) 

Heat  of  desorption  from  soil  to  air 
Mass  of  contaminant  trapped 
Duration  of  sampling 
Air-filled  porosity  (m3/m3) 

Water-filled  porosity 
Total  porosity  (m3/m3) 

Bulk  density  of  the  soil  (kg/m3) 

Kinematic  viscosity  of  air 
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